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Abstract 

We consider the many body quantum dynamics of systems of bosons interacting 
through a two-body potential N^^~^V{N^x), scaling with the number of particles 
N. For 0 < /3 < 1, we obtain a norm-approximation of the evolution of an appro¬ 
priate class of data on the Fock space. To this end, we need to correct the evolution 
of the condensate described by the one-particle nonlinear Schrddinger equation by 
means of a fluctuation dynamics, governed by a quadratic generator. 


1 Introduction 

In the last years important progress has been achieved in the mathematical understand¬ 
ing of the time-evolution of many body quantum systems. Here, we are going to consider 
the dynamics of systems of bosons, characterized by permutation symmetric wave func¬ 
tions. 

A bosonic system of N particles moving in three space dimensions can be described 
on the Hilbert space 

= {Uv e : HV’Aflb = 1 and . ..,a:^(Ar)) = -ipNixi, ■ ■ ■ ,xn)} 

The evolution of an initial iI^n ^ is governed by the Schrddinger equation 

idttpN,t = H]siil:N,t ( 1 - 1 ) 

where, on the r.h.s., 77jv is the Hamilton operator of the system. Restricting our attention 
to two-body interactions, the Hamilton operator takes the form 

N N 

Hn = + Vext{xj)) + X'^Vixi- Xj) (1.2) 

j=l i<j 
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where 14xt is an external potential, V is an interaction and A € M is a coupling constant. 
Notice that the unique global solution of (1.1) is given by = exp(— iffArt)'0Ar,o- 
For systems of interest in physics and chemistry, the number of particles involved in 
the dynamics is huge, ranging from N ~ 10^ up to values of the order N ~ 10^^. For 
this reason, despite the fact that (1.1) is a linear equation, it is typically very difficult 
to extract useful information. One of the main goals of non-equilibrium statistical me¬ 
chanics is therefore the derivation of effective equations approximating the solution of 
(1.1) in the interesting regimes. 

The simplest non-trivial limit, in which it is possible to obtain an effective approx¬ 
imation of (1.1) is the mean field regime, where A/ S> 1, |A| ^ 1 and N\ remains fix, 
of order one. These conditions guarantee that particles interact through a large number 
of weak collisions, whose total effect is comparable with their inertia. To investigate 
the time-evolution in the mean held regime, we set A = 1/N in (1.2). We obtain the 
Hamiltonian 

N N 

Hf = j;(-A,, + Hext(x,)) + - V{xi - Xj) (1.3) 

j=l i<j 

and we study the corresponding evolution '4)N,t = exp(—zFf™^t)V' 7 v Let us assume that 
the initial data is approximately factorized, i.e. that for a E L^(M^). 

Because of the mean-held nature of the interaction, we can expect that, for large A, 
factorization is approximately preserved by the time-evolution. In other words, we can 
expect that 'ipN,t — for a new (ft E L^(M^). Under this assumption, it is simple to 
show that must satisfy the self-consistent Hartree equation 

idtift = (-A -h Ve^t)v’t + iV * \ipt\‘^)(pt (1.4) 

where the cubic nonlinearity rehects the two-body interactions. 

To obtain a precise statement about the convergence of the many-body evolution 
towards the Hartree dynamics, we introduce the notion of reduced densities. The one- 
particle reduced density associated with the solution of the Schrodinger equation is 
dehned as the non-negative trace-class operator on L^(M^), with the integral kernel 

y) = j dX2... dXN X 2 ,..., XAr)'0Ar,t(y, X2, ■ ■ ■ , Xjv) 

normalized so that try^^^ = 1 ( 7 ^ ^ is obtained by taking the partial trace of the orthog¬ 
onal projection |'0Ar,t)(V’Af,t| over the degrees of freedom of the last (N — 1) particles). 

(k) 

Similarly, we can also introduce the A:-particle reduced density 7 jy ^ associated with 
for every k = 2,3,... , N. 

It turns out that the reduced densities 7 )^ ^ provide the appropriate language to dis¬ 
cuss convergence of the many body evolution towards the Hartree dynamics in the mean 
held regime. In fact, under appropriate conditions on the interaction potential V (in¬ 
cluding the case V{x) = ±|x|“^ of a Coulomb interaction), one can show that, for every 
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family of initial data V’A'',o G with |<y7)(y^| (approximate factorization at 

time t = 0), we will have 

'yNl^ ( 1 - 5 ) 

as N ^ oo, for all hxed t G M. Here (ft is the solution of the Hartree equation (1.4), 
with the initial data ^pt=a = ^- In fact (1.5) can be extended to get convergence of the 
fc-particle reduced density, for any hxed A: G N. The hrst results in the direction of (1.5) 
have been obtained in [28, 23, 38]. More recently, much work went into the proof of (1.5) 
in the case of singular interaction potentials; see [21, 4, 15, 22, 36, 1, 29, 30, 3, 11, 12, 9, 2]. 

After identifying the limiting effective dynamics (the one governed by the Hartree 
equation (1.4), in the mean held regime), it is natural to consider huctuations around it. 
To study huctuations, it is very useful to switch to a representation of the many particle 
system on the bosonic Fock space 

J- = 0L2(R3^,dxi...dx„) 

n>0 

On T, we can describe states with a variable number of particles. The vector 'k = 
{'00) V’l! • • ■ } £ ^ describes a state having n particles with probability H^nlli: 
n G N. For / G L^(]R^), we let a* if) and a(/) denote the usual creation and annihilation 
operators acting on T. We also introduce operator-valued distributions a*, ax creating 
and, respectively, annihilating a particle at x. They satisfy the canonical commutation 
relations 

[a^, 0*1 = d{x - y), [ax, ay] = [o*, a*] = 0 . 

In terms of these operator-valued distributions, we dehne the Hamilton operator on 

by 

= J dxa*x{-Ax+ Ve^tix))ax +j dxdyV{x - y)a*xa*yayax (1.6) 
Since commutes with the number of particles operator 

Af = J dx a^ax 

the corresponding time evolution preserves the number of particles. In particular, if we 
choose an initial data of the form T = {0,... , 0, ^iv, 0,... } with exactly N particles, its 
evolution will coincide precisely with the one generated by (1.3). 

The advantage of working in the Fock space, rather than in the A^-particle space 
is the freedom in the choice of the initial data. We are interested in the 
evolution of coherent initial data. For / G L^(M^), the coherent state with orbital / is 
given by W{f)Q G where H = {1,0, 0,... } is the vacuum state (with no particles) 
and 

VF(/)=exp(a*(/)-a(/)) 

is the Weyl operator with orbital /. Simple computations show that 

W(/)ll = e-ll^ll=/2|i,/,^,...| (1.7) 
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and that the expected number of particles is given by 

{w{m,Mw{m) = \\f\\i (1.8) 

Motivated by (1.8), we study the many body evolution generated by (1.6) for ini¬ 
tial coherent states of the form W{y/Nip)Q, with tp S L^(R^) such that ||v3||2 = 1 (this 
normalization guarantees that the expected number of particles in is equal 

to N). Since factorization is believed to be approximately preserved by the mean field 
dynamics, we expect that the evolution of the coherent state W{VNif)^l can be ap¬ 
proximated by the evolved coherent state W{y/Nipt)Q, where pt is the solution of the 
Hartree equation (1.4). 

We define the fluctuation dynamics 

= W{^/Niptye~^^^^W{^/Np) (1.9) 

and we set Then, we have 

e-^'^N'^WiVNip)n = W{VN(pt)(t ( 1 . 10 ) 

Hence, to prove that the full evolution of the initial coherent state (the l.h.s. of the last 
equation) is approximately coherent, we need to show that is close to the vacuum (in 
fact, it is enough to prove that the expected number of particles in is much smaller than 
N, since the evolved Weyl operator W{y/Nipt) creates a condensate with approximately 
N particles in the state pt)- To this end, it is useful to observe that 

z W(t) = ctmfit) 

with the generator 

T^^(i) = J iixa*(-Aa; Vext{x))ax + j dx{V * \(pt\‘^){x)alax 
+ J dxdyV{x - y)pt{x)Tpyy)alay 

+ ^ j dxdyV (x - y) [pt{x)pt{y)ala*y + h.c.] (1.11) 

+ j dxdyV{x - y)al [pt{y)a*y + h.c.] 

+ ^ / dxdyV(x - y)a*ayayaa, 

Notice that the terms on the third and fourth line do not commute with the number 
of particles operator J\f. This implies that the fluctuation dynamics does not 

preserve the number of particles (this is of course no surprise; the number of excitations 
of the condensate is expected to increase during the dynamics). Nevertheless it turns 
out that the expectation of N cannot increase too fast. In fact, using the expression 
(1.11), one can show that 

{iuMit) = < Ce^\^\ ( 1 . 12 ) 
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uniformly in N. This estimate on the growth of the expectation of M in the state can 
be used to show that ( 1 . 10 ) remains close to a coherent state, in the sense of the reduced 
densities. If we denote by 7 ^^^ the reduced density of the full evolution of W{-\/Nip)^ 
(more generally, of a state of the form W{y/Nip)^, with ^ having only few particles), 
( 1 . 12 ) allows us to show that 


tr 





< 


N 


(1.13) 


and that similar estimates hold for the fe-particles reduced density, for all A: G N. The 
study of the dynamics of coherent states has been initiated in [28, 23]. More recently, it 
has been further developed in [36, 11], leading to a proof of (1.13). 

A part from bounds like (1.13) on the rate of convergence of the reduced densities, 
this approach also allows us to study the fluctuations around the Hartree dynamics, in 
the limit of large N. From (1.11), we expect that, for N 00 , the fluctuation dynamics 
converges towards a limiting dynamics defined by 


with the generator 


0 <)= / 


= / dxali—A. 


IClxClx 


+ 

+ 


+ 14xt(a;))ax + J dx{V * \^pt\^){x)c 

J dxdyV{x - y)ipt{x)Tp^{y)alay (1.14) 

^ j dxdyV{x-y) [ipt{x)(pt{y)alal + h.c] 


independent of N. In fact, under appropriate assumptions on the interaction potential, 
one can indeed show that 




mf/ 


Vn 


(1.15) 


By definition, this implies that 


\e-i'H^'twiVN(p)n - wiVN(pt)uzHt)n\\ < 


Ce^|i| 

Vn 


(1.16) 


Hence, taking into account the limiting fluctuation dynamics we provide a norm 

approximation to the full many body evolution (hence, a stronger approximation com¬ 
pared with the one furnished by the evolved coherent state W{y/N(pt)^, which is only 
valid in the sense of the reduced densities). The convergence (1.15) has already been 
observed in [28] for smooth interactions and then in [23] for a larger class of potentials. 
More recently, it has been established, in a slightly different form, in [26, 27, 8 ]. Using 
(1.15) and bounds like (1.12) on the growth of the expectation of the number of particles 
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(and of its higher moments) with respect to the evolution one can prove a central 
limit theorem for sums of one-particle observables evolved through the full interacting 
many-body dynamics; see [5, 7]. 

Instead of considering fluctuations around the Hartree dynamics for coherent initial 
states on the Fock space, it is possible to analyze directly the mean field evolution in 
the A^-particle Hilbert-space for approximately factorized initial data. To this 

end, it is convenient to introduce the time-dependent map 


(1.17) 

where J-t denotes the bosonic Fock space, constructed over the orthogonal complement 
in L^(R^) of the one-dimensional space spanned by the solution ipt of (1.4) and = 

,..., , 0,0,... } if 

^ -t- (g)^ -^ ^s‘f + 


where (gs denotes the symmetric tensor product. When applied to the many body evo¬ 
lution of an approximately factorized initial data ipNfi, the isometric 

map UN,t eliminates the particles in the condensate tpt and let us focus on the fluctua¬ 
tions. It has been shown in [31], inspired by ideas developed in the time-independent 
setting in [32], that there exists a Fock space unitary evolution s) with a quadratic 

generator T[J^(t) such that 


UN,te Q)uNfi'lpN 


mf/ 


< CiV-i/2gX|t| 


(1.18) 


Notice that is similar but does not coincide with the limiting generator (1.14) (the 
difference between the two generators is due to the requirement, in the definition of UM,ti 
that fluctuations are orthogonal to (/?t). 

A more subtle and physically interesting regime, in which it is possible to approximate 
the many-body evolution by an effective dynamics, is the Gross-Pitaevskii regime. In 
the Fock space representation, the Hamilton operator is given by 

= J dxal{-A^ + Vextix))ax + ^ j dxdy N'^V{N {x - y))a*a*ayaoc (1.19) 

where P > 0 is a smooth, short range potential. It turns out that, in this case, the 
many-body Schrodinger evolution can be approximated by the time-dependent Gross- 
Pitaevskii equation 

idt^GP,t = (—A -I- Pext)</?GP,t + Sttoo |‘/^GP,tP<7’GP,t (1.20) 


where oq is the scattering length of the potential V. Recall that the scattering length is 
defined through the solution / of the zero-energy scattering equation 


-A + -P 


/ = 0 


( 1 . 21 ) 
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with the boundary condition f{x) 1, for |a:| —)• oo. For x outside the support of V, 
we have 

fix) = 1 - M 

\x\ 

where the constant oq is defined to be the scattering length of V. Equivalently, we can 
dehne oq through the integral 


Sttoo 


V{x)f{x)dx 


( 1 . 22 ) 


From the point of view of physics, the Gross-Pitaevskii regime is very different from 
the mean field limit since here particles interact rarely (only when they are very close, 
at distances of the order but when they do, the collisions are very strong. As a 

consequence of the strong interactions among the particles, the many body wave function 
develops a short range correlation structure which can be described by the solution / 
of the zero-energy scattering equation (1.21) and is responsible for the emergence of the 
scattering length in (1.20). 

A first derivation of (1-20) starting from many-body quantum mechanics has been 
given in [17, 18, 19]. Later, an alternative approach has been proposed in [35]. Bounds 
on the rate of the convergence towards the Gross-Pitaevskii dynamics have been then 
obtained in [6], making use of an appropriately modified version of the coherent states 
method illustrated above. The main problem one has to face to apply the coherent 
states approach to the Gross-Pitaevskii regime is the formation of correlations among 
the particles, which cannot be described by coherent states (the development of a short 
scale correlation structure in the Gross-Pitaevskii regime has been studied in [14, 10]). 
To circumvent this problem, one has to modify the coherent states through appropriate 
Bogoliubov transformations. Following [6], we define the function 

k?^{x-,y) = -N{1 - f{N{x - y)))(PGP,t(.x)cp^p t{y) 

where / is the solution of (1.21) and where T’qp ^ is the solution of a modified, N- 
dependent, version of (1.20), with the local nonlinearity replaced by an Hartree term, 
given by convolution with N^V{N.)f{N.) (since by (1.22), N^V{N.)f{N.) SttuqS, it 
is easy to bound the difference between ^ and (/Jgp,*)- Using we construct the 
unitary operator 


Ugp,* = exp 


- / dxdy[kf^{x\y)a. 


* 0 * 

Uy 



on T". Tgp ,4 is a Bogoliubov transformation; it acts on creation and annihilation opera¬ 
tors by 


UGP.ta(/)UGP,t = a(cosh^GP /) -h o*(sinh^GP /) 
?GP,t«*(/)^GP,t = a*(cosh^GP /) -h a(sinh^Gp7) 
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The operator Tqp,* can be used to implement the short scale correlation structure char¬ 
acterizing the solution of the many-body Schrodinger equation in the Gross-Pitaevskii 
limit. 

We consider the initial data W(\/]V<^)Tgp,o ^ (more generally, we can consider states 
of the form W{y/N (/?)Tgp,o?j with ^ containing only a bounded number of particles). We 
expect that the many-body evolution of such an initial data still has the same form. To 
confirm this fact, we define the fluctuation vector by requiring that 

e-^^^^W{VNip)TGP,on = ,)rGP,t 6 (l-23) 

From (1.23), to show convergence of the many-body dynamics towards (1.20) it is 
enough to prove that the fluctuation vector defined by (1.23) remains close to the 
vacuum. Since with 

= T^P,tW*{^/Nip^p^t)e-^'^^^W{VNip)TGP,o 


the problem reduces to show a bound for the growth of the number of particles with 
respect to the fluctuation dynamics . 

Such a bound has been established in [ 6 ], making use of certain cancellations in the 
generator of produced by the introduction of the Bogoliubov transformation TGP,t- 
As a consequence, it was proven in [ 6 ] that the reduced density 7 ^^^ of the full evolution 
of the initial data W{Vn ip)TGP, 0 ^ (or, more generally, of initial data having the form 
W{y/Nip)TGPfi^, for ^ G W with a bounded expectation for Af and T-Ln) satisfies the 
bound 


tr 




<CN exp(ci exp(c 2 |f|)) 


(1.24) 


where <pt is the solution of the Gross-Pitaevskii equation (1.20). 

In the mean held regime that we discussed above, the coherent states approach could 
also be used to describe huctuations around the limiting equation. In particular, it al¬ 
lowed us to identify a limiting huctuation dynamics with a quadratic generator and to 
apply it to obtain a norm bound of the form (1.16) for the many-body evolution. After 
establishing the estimate (1.24) for the rate of convergence of the one-particle reduced 
density, it is therefore natural to ask whether we can use the same approach to describe 
huctuations around the Gross-Pitaevskii equation in the limit of large N. Unfortunately, 
it turns out that the in the Gross-Pitaevskii regime, one cannot approximate the huctua¬ 
tion dynamics by a quadratic evolution in norm. Although one can control its effect 
on the growth of the number of particles (needed to prove (1.24)), the cubic and quartic 
components of the generator of (cubic and quartic in the creation and annihilation 
operators) are not negligible in the limit of large N. 

Instead of considering huctuations of the time-evolution around the time-dependent 
Gross-Pitaevskii equation, it is also possible to approach this problem from a static, 
time-independent, point of view. To this end, one can trap the system in a hnite volume 
(either by imposing boundary conditions or by turning on an external potential) and one 
can study the difference between the many-body ground state energy and the minimum 




of the Gross-Pitaevskii energy functional or, more generally, the energy of low lying 
excitations (the fact that Gross-Pitaevskii theory describes, in leading order, the ground 
state properties of the many body system has been established in [34, 33]). In the 
mean field setting, this program has been carried out in [37, 24, 32, 13], where it was 
proven that the excitation spectrum is determined by a quadratic Hamiltonian similar to 
(1.14). This suggests that, in the mean field regime, a good approximation for the many 
body ground state has the form W{y/Nip)TVl^ where minimizes the Hartree energy 
functional and where T is a Bogoliubov transformation (the exponential of a quadratic 
expression in creation and annihilation operators), needed to diagonalize the quadratic 
Hamiltonian. Similarly, good approximation for low-lying excited states have the form 
Wa*{pg\)... a*{gk)^, for appropriate A; G N, A^o = ~ ^ and orbital gi,... ,gk 

orthogonal to p (in fact, to produce states with a fixed number of particles, it is better 
to work with a map defined similarly to (1.17), rather than with the Weyl operator 
W{^/Np)-, see [32] for details). 

Although the excitation spectrum in the Gross-Pitaevskii regime should still be close 
to the spectrum of a quadratic Hamiltonian, a good approximation for the ground state 
cannot have the form W{-\/Np)TQ. for a Bogoliubov transformation T (analogously, 
excited states cannot be approximated by vectors like W{y/Np)Ta*{gi)... a*{gu)^)■ 
This follows from [20], where it has been shown (in fact, in a more general setting) 
that the minimum of the energy over all states of the form W{y/Np)TQ. (with T being 
the exponential of a quadratic expression) remains strictly above the true ground state 
energy, with an error of order one (an upper bound to the correct ground state energy, 
up to an error that, in the Gross-Pitaevskii regime, vanishes in the limit of large N, has 
been obtained in [39]; this result is consistent with the Lee-Huang-Yang prediction). 

In this paper, we are going to consider an intermediate regime, lying between the 
mean field and the Gross-Pitaevskii limits. For 0 < /3 < 1, we define the Hamilton 
operator 

Tfw = y dxV^alV^a^ + ^ / ^xdy N^^V{N^{x - y))<a;aj,a, (1.25) 

acting on the bosonic Fock space F. To simplify a bit the computations, we neglect 
here the external potential (but it would be easy to modify our analysis to include one). 
The Hamiltonian (1.25) can be thought of as an interpolation between the mean field 
Hamiltonian (1.6), obtained with /3 = 0, and the Gross-Pitaevskii Hamiltonian (1.19), 
recovered with (3 = 1. 

For 0 < /3 < 1, the many body evolution develops weaker correlations, compared with 
the Gross-Pitaevskii regime. As a consequence, on the level of the reduced densities, the 
many body evolution generated by (1.25) can be approximated, in the limit N ^ oo 
and for all 0 < /3 < 1, by the nonlinear Schrodinger equation 

idtpt = -A<y9t bolptl'^pt (1-26) 

with bo = f V{x)dx (notice that Svrao < &o for all short range and repulsive V). 

While (1.26) is enough if we are interested in the limiting behavior of the reduced 
densities, to study fluctuations and to obtain a norm approximation for the many body 
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evolution we need a more precise ansatz, taking into account the (weak) two-body cor¬ 
relations. Instead of working with the solution of the zero-energy scattering equation 
(1.21), we find more convenient here to fix £ > 0 and to consider the ground state /jv,£ 
of the Neumann problem associated with the potential .) on the sphere of 

radius i centered at the origin. In other words, we choose /tv,^ as the solution of the 
eigenvalue problem (2.4) associated with the smallest possible eigenvalue normal¬ 
ized so that fN,e{x) = 1 for |a;| = i and continued to by requiring that fN,iix) = 1 
for all \x\ > i. We use to describe correlations among particles in the condensate. 
Accordingly, we consider the A^-dependent Hartree equation 

P)¥^f . (1.27) 

As N —^ oo and for all 0 < /3 < 1, approaches the solution of the nonlinear equa¬ 
tion (1.26). We will see, however, that (1.27) furnishes a better approximation for the 
dynamics of the condensate wave function, because, through the factor Jn/, h takes 
into account the correlations among the particles (which, despite being weak, are not 
negligible in the analysis of the fluctuations). 

Our goal is to study the fluctuations around (1.27), to prove that their dynamics 
has a quadratic generator in the limit of large N, and to use the limiting fluctuation 
dynamics (with the quadratic generator) to obtain a norm approximation of the many 
body evolution generated by (1.25). 

First of all, we need to take care of the correlation structure. We proceed similarly as 
in [6], introducing a family of Bogoliubov transformations. Using the Neumann ground 
state Jn/, we set = 1 — In/ (so that ujN,t{x) = 0 for all |xl > i) and we define 

y) = -NuN/ix - y){(pf{{x + y)/2)f (1.28) 

where is the solution of (1.27). It turns out (see Lemma 2.1) that ||A;Ar,t ||2 is bounded, 
uniformly in N, and therefore (1.28) is the integral kernel of a Hilbert-Schmidt operator 
that we denote again by few,*- Using we dehne the Bogoliubov transformation 


TN,t = exp 


J dxdy kisf^tix',y)a*a* — h.c. 


(1.29) 


We consider initial data of the form W, with a “close” to the vacuum 
U (in the sense that the expectation and the variance of the number of particles and of 
the kinetic energy operator in the state can be bounded uniformly in N). We write 
the evolution of such initial data as 




where 

iN,t = 7/7v(i; 0)^N 

with the fluctuation dynamics 

UN{t-,0) = T;f^tW*iVNp^)e-^^^^W{VNip)TN,o 
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Our goal is to approximate 0) by an evolution with a quadratic generator. 

To this end, we define the phase 

r]Nit) = N j dxdyN'^^V{N^{x-y)){l/2-fN/{x-y))\^t {x)\^\^Pt {y)\^ 

+ j dxdy\V^ smhk^^^{x, y)\‘^ + j dx{N^^V{N^.) * \‘^){x){s^, s^) 

+ J dxdyN^Pv{NP{x-y))^^{x)(p^{y){s^,s^) 

+ Re I dxdyN^^V{N^{x-y))ip^{x)^^{y){s^,c^) 

+ ^/ dxdy {N^ {x - y)) [|(sf, c^)|^ + |(s^, + (s^, s^)(sf, sf)] 

(1.30) 

and the time-dependent and iV-dependent quadratic generator 

C2At) = + 45(t) 

N f 

+ Y dxdyuN/ix-y) 

X ((a^ + y)/2)A(y9f ((x y)/2) + \V(p^{{x + y)/2)|2)o>* -h h.c.] 

+ NXi^n f dxdy l(|x - y\ < i) [(v^f ((a; + v)/‘^)falal + h.c.] 

(1.31) 

with 

= jdx {N^PV{nP.) * W^?){x) 

X [a*(cf)a(cf) + a*(sf)a(sf) + a*(cf)a*(sf) + a(sf)a(cf)] 

J dxdyN^^V{N^{x-y))(p^{x)(f^{y) 

X [a*{c^)a{c^)+a*{s^)a{s^)+a*{c^)a*{s^)+a{s^)a{c^)] 

+ ^J dxdyN^^V{N^{x-y))(pf^{x)(p^{y) 

X [a* (cf )a(s^) a* {c^)a{s ^) + a(sf )a{Sy )] 

+ ^J dxdyN^PV{Nd{x-y))(p^{x)^^{y) 

X [a* {s^)a{c^) + a* {s^)a{c^) + a* (sf )a* (s^)] 

+ ^ j dxdyN^^V{Nl^{x-y))ipf{x)(p^{y) [a*{p^)al + a*{c^)a*{p^)] 

+ 11 dxdyN^f^V{Nf^{x-y))ipf{x)ipf{y) [a{p^)ay + a{c^)a{p^)] 

(1.32) 
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and 



(ixVa;0* Vxflx + I dx 




ala{-A^p^) +a*(-A3;p")aa; + a*(Va;p")a(Va;P^'') 


N\ 


„N\ 


^N\ 


+ V,a*ik,)VMkx) + a*{-A,r^)a{k,) + a*{s^)a{-A,r^) 
+ a*{-A^p^)a*{kx) + a{kx)a{-Ap^) + ala*{-A^r^) 

+ a{-A^r^ )a^ + a*{p^ )a* (-) + o(-A^r^ )a(pf) 


(1.33) 


Here we used the shorthand notation (y) = cosh^j^ ^(y, x), {y) = sinh^^^ ^(y, x), 

Px (y) = (coshfcj^_^ -l){y,x) and (y) = (sinh^j^^ -kN,t){y,x). 

We denote by V( 2 ,N{t]s) the time evolution generated by C 2 ,n- In other words, 
^ 2 ,Af(i; s) is the two-parameter group of unitary operators solving the Schrodinger equa¬ 
tion 

idtU 2 ,N{t] s) = C 2 ,N{t)ld 2 ,N{t] s) (1-34) 

with the initial condition hl 2 ,N{s\ s) = 1 for all s G M. The existence of such an evolution 
U 2 ^n can be established as in [23] (notice here that all kernels entering the generator 
C- 2 ,N{t) are smooth for every finite A; they only develop singularities as A —)■ oo). 

In the next theorem, our first main result, we show that, when acting on appro¬ 
priate vectors ^Af(i;0) can be approximated by Z^ 2 ,Ar(t; 0 ), up to the (physically 
uninteresting) phase rjN- 


Theorem 1.1. Let V >0 be smooth, spherical symmetric and compactly supported. Fix 
0 < /3 < 1 and let a = min(/3/2, (1 — /3)/2). Let ip G A^(M^). Fix £ > 0. Consider a 
sequence G such that ||^Af|| = 1 and 

{Cn, [Af‘^ + JC"^ + Vat] ^n) < C (1.35) 


uniformly in A. FIere we introduced the notation 

IC = I dxV,<V,a,, Vn = ^ J dxdyN^^V{N^{x - y))ala*yaya,, (1.36) 

for the kinetic and the potential energy operators entering in LLn = + Vn- Then there 

are C,ci,C 2 >0 such that 

< CN~°‘ exp(ci exp(c 2 |t|)) 


for all t gM. and all A large enough. 

Notice that a result similar to Theorem 1.1 has been recently obtained in [25], however 
only for 0 < /? < 1/3. The main difference with respect to [25] is the fact that here we 
already introduce the Bogoliubov transform Tat^o at time t = 0; this allows us to cover 
all /3 < 1. 
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Theorem 1.1 shows that, in the limit of large N, the dynamics of the fluctuations 
around the nonlinear Schrodinger equation (1.27) can be approximated by the evolution 
having the quadratic generator (1.31). It is also possible to approximate the 
dynamics of the fluctuations by a limiting evolution, again with a quadratic generator, 
but now independent of N. To this end, we start by noticing that, in the limit of large 
N, the integral kernel (1.28) approaches the limit 

kt{x] y) = + y)/2) (1.37) 


where ypt is the solution of the independent nonlinear Schrodinger 
while 


asymp 


(x) 


f bo 


-f 

1 

< Stt 

Jx| 

1 0 




for |x| < i 
otherwise 


equation (1.26) 
(1.38) 


is the pointwise limit of Nujn^i{x) for N oo. Like /cAr.tj also kt S 
the integral kernel of a Hilbert-Schmidt operator on L^(M^). It defines 
transformation 


Tt = exp 


^ j dxdy ktix-y)a*a*y-h..c. 


L^(M^ X M^) is 
the Bogoliubov 

(1.39) 


With this notation, we can formally take the limit —)■ oo in (1.31) and define the 

time-dependent quadratic generator 


C2,ooit) = {idtTm + Cf\t) + 6P{t) 

+ ^J dxdy (x - y) 

X { [(pt{(x + y)/2)A(pt{{x + y)/2) + \V(pt{{x + y)/2)|^] -h h.c.) 

+ J dxdyl{\x -y\<i) + y)/2)a>* -h h.c.] 

(1.40) 

where and are defined as and but with and k]\f,t replaced by 

(ft and kt (and thus with replaced by Cx, Sx,Px,rx, dehned analogously). 

We denote by Vl 2 ,oo{t]s) the two-parameter unitary group generated by (1.40). In 
other words, we define 7 / 2,00 through the Schrodinger equation 


Z(9t/72,oo(/) s) — 712,00(/)7/2, 00 (/) s) 


with the initial condition U 2 ,oo{s‘,s) = 1. Also here, the existence of such an evolution 
7 / 2,00 is guaranteed by the work [23] (the only singularity entering £ 2,00 is the Coulomb- 
type singularity of 

Theorem 1.2. Let V >0 be smooth, spherical symmetric and compactly supported. Fix 
0 < /3 < 1 and set a = min(/3/2, (1 — /3)/2). Let p £ Fix £ > 0. Consider a 

sequence such that H^Afll = 1 CLnd 

{Cn, + IC'^ + Vat] Cat) < C 
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uniformly in N. Then there exist C, ci, C 2 > 0 sueh that 

\e-^'^^^W{VN^)TN,o(N - e-^foVN(s)ds W{VNip^)TN,tU2,oomNf 

< CN~‘^ exp(ci exp(c 2 |t|)) 

for all t and all N large enough. 

Theorem 1.2 shows that, for all 0 < /3 < 1, the dynamics of the fluctuations around 
the nonlinear Schrodinger evolution (1.27) is governed by the quadratic generator (1.40). 


Aeknowledgements. The authors would like to acknowledge support by the Swiss Na¬ 
tional Science Foundation through the SNF Project “Effective equations from quantum 
dynamics”. S. Cenatiempo acknowledges the support of MIUR through the FIR grant 
2013 “Condensed Matter in Mathematical Physics” (code RBFR13WAET). C. Boccato 
and B. Schlein also gratefully acknowledge support by the CRC-1060 “The Mathematics 
of emergent effects”. 


2 The Fluctuation Dynamics 

We work on the bosonic Eock space 

^ = C©0l2(m3’") 

n>l 

where L^(M^”) denotes the space of all ijjn S T^(M^"') such that 

'0n(^7rl) • • • ) X-Kn) — '0n(^l) • • • j ®n) 

for all permutations tt £ Sn- We use the notation = {1,0,0 ... } for the vacuum vector 
in T. 

Eor g G L^(M^), we introduce the creation operator a*{g) and its adjoint, the anni¬ 
hilation operator a{g), by 

1 

{a*{g)'^)^^\xi, ...,Xn) = g{xj)'ip^'^~^\xi ,.. .,Xj-i,Xj+i, ...,Xn) 

i=i 

(a(g)T)("')(xi, ...,Xn) = Vn + l j dx g{x) il}^'^^^\x,xi,... ,Xn) 

for all T = G Creation and annihilation satisfy canonical commutation 

FclcttlODG 

[«(/),«* (s')] = (/, S'), [a(/), a(s')] = [a* if), a* (S')] = 0 

We will also make use of the operator-valued distributions a*, , defined so that 

«*(s) = j dxg{x)a*, aig) = j dxg{x)ar, 
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Although creation and annihilation operators are unbounded operators, they can be 
bounded with respect to the number of particles operator, defined by 

(AA^)W = 

for all 'h = G F or, equivalently, in terms of the operator valued distributions 

by 

Af = J dx a*ax ■ 

Using this last expression, it is easy to check that 

I|o'(9)'I’II<M|2|IW + i)'''"'I'I| 

On F, we are interested in the time-evolution generated by the Hamilton operator 

'Hn = J dxVxalVxa^ + ^ N^^V{nI^{x - (2.1) 

for 0 < /3 < 1 and for a smooth potential U > 0 with spherical symmetry and with 
compact support. Notice that commutes with the number of particles operator AA; 
hence, the time evolution preserves the number of particles. On the n-particle sector, 
(2.1) acts as 

n 1 ^ 

= Y, -A,^, + - y)) 

j=l i<j 

We are going to study the time-evolution generated by (2.1) on coherent states. For 
g G L2(M3), we define the Weyl operator 

W{g) = exp(a*( 5 ) - a{g)) = e-<s) ( 2 . 2 ) 

Weyl operators are unitary W*{g) = W{—g) = W~^{g) and they act on creation and 
annihilation by shifts, i.e. for any f,g & L^(M^), we have 

W{gra{f)W{g) = a{f) + {f,g), 

Wigra*if)Wig) = a*if) + {gJ) ^ ^ 


For ip G L^(M^) with \\p >\\2 = 1, the coherent state W{s/N^ p)Q, describes a condensate 
with an average of N particles, all described by the orbital (recall (1.7). To obtain 
a norm approximation of the many body evolution, we have to implement the correct 
short scale correlation structure on top of W{-\/Np))Vl. To this end, we fix £ > 0 and we 
define /tv/ to be the ground state of the Neumann problem 


-A + -^N^^V{N^x) 


fN,e — ^N,efN,e 


(2.4) 


on the sphere jxl < i, normalized so that fN,i{x) = 1 for \x\ = ^. In the next lemma, 
whose proof is deferred to Appendix A, we collect some important properties of the 
solution of (2.4). 
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Lemma 2.1. Let V be smooth, positive, spherically symmetric and compactly supported 
with bo = f Vdx. Let /jv,£ be the ground state of the Neumann problem 

( — ^ + 2^^^ {^^'))fN = ^N,efN,i (2.5) 

on the sphere of radius i, with the boundary conditions 


fN,i{x) = 1 drfN/{x) = 0 

for all X G with |x| = i. For N sufficiently large (such that RN~^ <£) we have: 

i) 


^N1 — 


3bo 


8ttNP 


< 


C 




a) There is a constant 0 < cq < 1 such that, for all jxj < i, 

Co < fN,e{x) < 1 


( 2 . 6 ) 


(2.7) 


Hi) Let UN^i = 1 — fN,£- There exists a constant C > 0 such that, for all \x\ < i, 

- Af(|x| + JV-g) S + jv-2^) (2-8) 


We continue f]y,e to a function on by setting fN,e{x) = 1 for all lx| > £ and we 
define = 1 — fjy^^ (so that a;jv,r(a^) = 0 for all |x| > £). To generate correlations in 
the condensate W{'/N(f)Ll at time t = 0, we define 

kN,o{x; y) = -Nujn,(.{x - y)gP‘{{x + y)/2) (2.9) 


With Lemma 2.1 it is easy to check that H^Ar^olb < C, uniformly in N (the constant 
depends on £, though). Hence, (2.9) is the integral kernel of a Hilbert-Schmidt operator, 
which we will denote again by k]\f,o- We use kN,o to define the Bogoliubov transformation 


rAr,o = exp 


J dxdykN,oix] y)ala* - h.c. 


acting as a unitary operator on T. 

We are going to study the time evolution of initial data of the form W{y/N(p)T]\[^of,N, 
for a sequence G T" satisfying (1.35) (this assumption guarantees that our initial data 
are dominated by the Weyl operator W{y/Nip), creating a condensate with an average 
of N particles in the state ip, and by the Bogoliubov transformation T7 v,0j creating the 
short scale correlation structure). 

To construct an ansatz to approximate the full evolution of such an initial data, we 
proceed as follows. First of all, we let the condensate wave function (p evolve, according 
to the Hartree equation 


idtp^ = -Apf + {N^^V{N^.)fN,, * |7>f P)<^f (2.10) 
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with the initial data Some standard properties of this equation (including 

global well-posedness in the energy space and propagation of higher Sobolev 

regularity) are reviewed in Proposition B.l. 

On top of the evolved condensate, we add the short scale correlation structure, ap¬ 
propriately modified to take into account the variation of the condensate wave function. 
Similarly to (2.9), we define 

kN,t{x-, y) = -NujN,iix - y){ipf{{x + y)/2)f (2.11) 


Like at time t = 0, the function turns out to be the integral kernel of a Hilbert- 
Schmidt operator. We use it to define the family of Bogoliubov transformations 


TN,t = exp 


J dxdykN,t{x]y)ala*y 



( 2 . 12 ) 


Observe that the unitary operators act on creation and annihilation operators by 


Th,ta{9)TN,t = a(coshfc^ ^ g) + a* (sinh*,^ ^ g) 

TN,ta*ig)TN,t = a*{cosh.kj^y g) + o(sinhfc^_^ g) 

where 

coshfc^^ = 7^ ikN,tkN,tT , sinhfc^_^ = Y (2n + If! ‘ 

We define the fluctuation vector at time t, requiring that 

e-^^^^W{^/Nip)TN,o^N = W{VN^f)TN4N,t (2.14) 

Equivalently, we have = ^Af(L0)?Af with the fluctuation dynamics 

UN{t; s) = )Tjv,. (2.15) 

The fluctuation vector measures the distance from the modified coherent state 
W Following the analysis of [6], one can prove a bound of the form 

< (7exp(ci exp(c 2 |t|))(^Ar, [Af + /N + nN] ^n) 

for the growth of the expectation of the number of particles with respect to the fluctu¬ 
ation dynamics. As explained in the introduction, this bound provides an estimate of 
the form (1.24) for the convergence of the reduced density associated with the full many 
body evolution (2.14). 

Here, we want to go one step further, and prove a norm-approximation for (2.14), 
using (2.10) for the evolution of the condensate, (2.12) for the description of the corre¬ 
lation structure, and a unitary evolution with a quadratic generator to approximate the 
fluctuation dynamics. To this end, we will use the following proposition. 
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Proposition 2.2. Let P > 0 6e smooth, spherical symmetric and compactly supported. 
Fix 0 < /3 < 1 and set a = min(/?/2, (1 — /3)/2). Let ip G Fix £ > 0. Consider 

a sequence G such that ||^Ar|| = 1 and 

{^N, + Vat] ^n) < C 


uniformly in N. LetUN he the fluctuation dynamics defined in (2.15), andU 2 ,N be the 
quadratic evolution defined in (1.34). Then there exist C,ci,C 2 > 0 such that 


h{N{i',0)fN — e '' U 2 ,N{t; 0 )f,N 


<CN " exp(ci exp(c 2 |t|)) 


for all t and all N large enough. 

Making use of this proposition, we can easily prove Theorem 1.1. 
Proof of Theorem 1.1. We notice that 


- e-'/o mis)ds )rw,tW2,iv(i; 0)^w 

W{^/Npf^)TN,t [uN{t]0) - e-^So^^i^)<i^U2,N{t]0)\ Cn 
ZYAr(t;0) - f,N < C'7V“"exp(ci exp(c 2 |t|)) 


where, in the last inequality, we used Proposition 2.2. □ 

To prove Proposition 2.2 we will compute the generator of the fluctuation dynamics 
Un and we will compare it with the quadratic generator (1.31) of U 2 ,n- This is the 
content of the next two sections. 


3 The Generator of the Fluctuation Dynamics 

From the dehnition of the fluctuation dynamics (2.15), we find 

idtUN{t] s) = CNit)V(N{t; s) 


with the generator 

CNit) = {im^fiTr,,t + T*^^, [{idtW*iVNp^)) W{VNp^) 


(3.1) 


The goal of this section is to prove the following theorem, describing the properties 
of the generator CN{t)- 
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Theorem 3.1. Let V > 0 be smooth, spherical symmetric and compactly supported. Fix 
0 < /3 < 1 and set a = min(/3/2, (1 — /3)/2). Let (p G Fix i > 0. Let Un be the 

fluctuation dynamics defined in (2.15). Then 


^N{t) = VN{t) + ^2,Af(i) + Vat + fAr(t) 

where the phase r]N{t) and the quadratic generator C 2 ,N{t) are given by (1.30) and (1.31), 
respectively, Vn is defined in (1.36), and where the error term FnA) satisfies 

±£N{t) < 5 Vn + CM^/N + Cse^^^^iM + 1) 

±[Af,£Nit)] < 6 Vn + CAf^ /N + Cse^^^^iAf + 1) (3.2) 

±£N{t) < + CAf^/N + Cse^^^^iAf + 1) 

and 


\{^Pu£Nm2)\ < C7V-“e^l'l [{^1, {}C+M+l)^Pi) + {1P2, (/C2 + (Ar + 1)2)^2)] (3.3) 

for all 1)1,1)2 £ iF. Moreover, the quadratic generator C2,N{t) satisfies the bounds 

±{C2,N{t) - JC) < Ce^\^\{Af + 1), {C2,N{t) - JCf < Ce^\^\{Af + 1)^ 

±[Af,C2,N{t)]<Ce^\^\{Af + l), ±[Af‘^,C2,N{t)] <Ce^\^\{Af + lf (3.4) 

±£ 2 ,Nit) < Ce^l*l(Ar+ 1), \£2,Nit)\^ < Ce^l'l(Ar+ 1)^ 

To prove Theorem 3.1, we compute the different terms on the r.h.s. of (3.1). From 
(2.3) we obtain 

4 

iidtW*iVNp^))wiVNp^) + w*iVNp^)nNW{VNp^) = Y^g^^\t) 

1=0 


with 


gj^)(i) = N j dxdyN^^ViN^ix - y))(l/2 - - y))\p^ix)\^\p^iy)\^ 

G^^\t) = Vn [a*i{N^^V{N^)ujN,e * 

+ai{N^^V{Ni^.)u}N,(: * |<^f ) 

JdxV,alV,a, + J dx (N^^ViN^.) * \pV^)V)a>x 

j dxdyN^I^ViN^ix - y))ip^ix)Tp^iy)alay 

i j dxdyN^^ViN^ix-y)) [ 99 ^(x)<y9f (y)a*a* + (x)^f (y)a,rS 

oPit) = ^ j dxN^^V{N<^ix - y))al {pf iy)a*y +Tp^iy)ay) a,, 

^ j dxdyN^^ViN^^ix - y))ala*yayax 


+ 

+ 
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From (3.1), we conclude that 


CN{t) = [idtT*^^,] TM,t + 

j=0 

Next, we analyze the contributions to the generator £Ar(t) separately. 

For j = 0, we have 


n,tSS\t)TN,t = N j dxdyAr^^F(iV^(x-y))(l/2-/^,Ka;-y))l¥^f(^)Pl¥’f(y)l' (3-5) 


For j = 1, we just write 


T*N,tSN\t)TN,t = a*{{N^f^V{Nf^.)coN,e * l</^f P)<^f) + h.c. 


Tnj 


(3.6) 


since at the end the main part of this term will be canceled. 

In the next four subsections, we study the terms for j = 2,3,4 and 

the term {idtT^ At the end, in subsection 3.E, we combine all results to prove 

Theorem 3.1. 


3.A Analysis ot 

We denote by 

KL — dx 

the kinetic energy operator. We will use the shorthand notation CN,t = cosh^^j, SN,t = 
sinhfc^^j, pN^t = CN,t - 1, rN,t = SN,t - kN,t and, for a fixed x G M^, (y) = CN,t{y]x), 

Sx (y) = snAV’ Px (y) = PN,t{y; x) and (y) = rN,t{y; x). 

The goal of this subsection is to prove the following proposition, which is a conse¬ 
quence of Proposition 3.3 and Proposition 3.5 below. 

Proposition 3.2. Under the assumptions of Theorem 3.1, we have 

Tf,^,QfTM,t= I dxdy|V,sinh,^_,(x,y)l2 + J dxiN^^V{N^) * \pf^\^)ix){s^, s^) 

+ 1 dxdyN^^V{N^{x-y))pf{x)p^iy){s^,s^) 

+ Re j dxdyN^^V{Nf^{x-y))p^{x)p^{y){s^,c^) 

+ j dxdy{-A^kN,t{x,y))ala*y + j dxdy {-A^kN,t{x,y))a^ay 

+ ^ I dxdyN‘^^V{Nf^{x-y)) {x)p^{y)ala*y + (p^{x)Cp^{y)axay\ 

+ 

(3.7) 


20 




(3.8) 


with and (1.33) and, respectively, in (1.32). Putting 




,{K) 


(V), 


-2,AT 


2 ,N\ 


we find 


±(4'^(t)-/C)<Ce^l*l(Ar + l), 

'Af,cf),it)] <Ce^l*l(Ar+l), 


{cfl^{t)-}Cf < C'e'^l*l(Ar+ 1)2 


Ar2^44(t)J <Ce^N(Ar + i)2 


±£g4t)<Ce^l*l(AA + l), 

To prove Proposition 3.2, we split 

T^,tGS\t)TN,t = n,tK-TN,t+n^,G^^’^\t)TN,t 

with 


(3.9) 


q{2 


'^\t) = J dx {N^^V{N^.) * P)(x) 


CL^dx 


+ J dxdyN^^V{N^{x-y))^^{x)ip^{y)alay 


+ 


dxdy N^f^V {x - y)) [yjf (x)(^f (y)a*a* + iff {x)ipf {y)aa:ay 


(3.10) 


We start by analyzing From (2.13), we find 

T(‘ICTt= j dxdy\Vr,smhktiy,x)\‘^ 

+ J dxdy{-ll^xkN,t{x,y))ala*y + J dxdy {-A,j;kN,tix,y))axay 
+ ci^]^{t) 

with ^ 2 ^ as given in (1.33). 

Proposition 3.3. Under the assumptions of Theorem 3.1, the operator C^j^{t) defined 
in (1.33) satisfies the bounds 

±(45J(0 - < Ce^'''(AA + 1), (45 - + 1)^ 

■AA, 45 j <Ce^l^l(AA + l), ±[Af^ci^^{t)]<Ce^\\Af + lf ( 3 . 11 ) 


±4'^(t)<C7e^l^l(AA + l), 


45(^) " ^ Ce^'*'(AA + l)' 


The proof of Proposition 3.3 is based on the estimates contained in the next lemma. 
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Lemma 3.4. Let ji,j 2 S L^(M^ x M^) and denote ji^x{z) = ji{z,x) fori = 1,2. Set 



Then, we have 

|(V’,AiV’)| < C'||ji||2||(AA + 

|(V',^2^)|<C||ji||2|U2||2||(AA+l)V2^f 


(3.12) 


for all if a J-. Moreover, 

A\A^ + A^A\<C\\h\\l{M + lf 

A2A2 + ^ 2^2 — C'||jl||||j2||2 {Af + 1)^ 

Furthermore, let 

A 3 = J dxa*{Vxkx)a{Vxkx) 

where we put k^(y) = k]\f^t{y, x), with kN,t as defined in (2.11), with initial data ipN,t=o = 
ip £ Then we have 

A 3 < Ce^l‘l(AA+1), Ai < (76^1*1 (AT + 1)2 (3.14) 


and also the bound 

A^As + A 3 A^ < {Af + if (3.15) 

for the time derivative 

A3 = J dxa*{Vxkf)a{Vxkf) + j dxa*{Vxkf)a{Vxkf) (3.16) 

Proof. We start with (3.12). To this end, we observe that 

|(V',AiV')|< j dx\\afji^x)'ip\\\\axfj\\ 

< J C^a:||jl,,i;||2||(AA+ 1)^/2 ^||||q^^|| < l)l/2^||2 

by Cauchy-Schwarz. The bound for A 2 follows similarly. Now, we show (3.13). We begin 
with 


{'ip,AlAiil>) = j dxdy{'4},alafji^x)a^{ji^y)ay'il}) 

J dxdy\\afji^x)ax'if\\\\afji,y)ay'ilj\\ 

J <^a:dy||jd,x||2||i2,y||2||axAA^^^'0||||ayAr^/2^|| 

r r 1 1/2 r i- 

J dxdy\\ji^x\\l\\ayAf^^‘^'iff J dxdy\\ji^y\\l\\axAf^^‘^'i)\('‘ 


< 

< 


< 


< 


- 11/2 


^ IUIII 2 
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Analogously 

= j dxdy{'ilj,J{ji^x)axala\ji^y)'ip) 

< J dxdy{'4!,J{ji^x)alaxa^iji,y)'ilj)+ J dx{'ilj,J{ji^x)a\ji,x)'tp) 

< J dxdy\\ayJ{ji^x)ij\\\\axa^iji,y)ij\\ + \\jif\\{Af + 


which implies, as above, that {ip, AiAlip) < 2|| ji II 2 II (A^ + 1)V’|P- The estimate for 
A 2 A 2 + A 2 A 2 follows analogously. 

We switch now to the term A 3 . We write it as 


A3 = / dxdydzVxkN,t{v,x)S/xkN,t{z,x)alay = 


J dy a*{uy)c 


with 


Uy (z) = UN,t{z,y) = I dxVxkN,t{y,x)VxkN,t{z,x) 


To bound the x M^) norm of U 7 v,tj we notice that 

hN,t\\l = / dzdy 1 / dx VxkN,t{y-, x)VxkN,t{z, x)| 

= J dzdydxidx2 V xikN,t{y-,xi)V xikN,t{z,xi)V x2kN,t{v,X2)V x2kN,t{z,X2) 


Since 


\VxkN,t{x\ z)\ < [Iv^f ((® + z)/ 2 )\^ + \Vip^{{x + z)/ 2 )\^] 


X — z 


we find 


\\uN,t\\l < j dxidx 2 dy [|(y9f ((xi + y)/2)p + ((xi + y)/2)|2] 

1(1x1 -y\< ^)l(|x2 -y\<i) 


/ 


dz-: 


< J dxidx2 


|xi - ?/p|x2 - yp J \xi - z\‘^\z - X2\^ 
l(|xi| < •^)1(|X2| < i) 


|xi|21x2P1xi - X2I 


Hence, (3.14) follows from (3.12) and (3.13). As for (3.15), we proceed similarly, re¬ 
marking that 

if'fe,.) = /dxV.j„,fe..)V.E„,,(=,x) and 

d2iy,z)= / dxVxkN,t{y,x)VxkN,tiz,x) 
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are such that ||uf ^||2 < for i = 1,2 (one has to observe that one can avoid to take 

the terms Vip^ in the L°°-norm). □ 


Proof of Proposition 3.3. Since 

||'SAr,4||2, ||P7V,t||2, IkTV.tIb, \\'^xPN,t\\2, \\^xPN,t\\2, ||||2 < 

are all finite, uniformly in N (see section C), we note that the terms in (1.33), with 
the exception of the kinetic energy operator /C, have the form of one of the operators 
Ai, Af, A 2 , A 3 analyzed in Lemma 3.4. From (3.12), (3.14) it follows immediately that 

With Cauchy-Schwarz, we can bound the square of — JC) by the sum of terms like 

AiA* and A*Ai, with i = 1, 2, 3. Lemma 3.4 implies therefore that 

Since the commutators of N with the terms on the r.h.s. of (1.33) are either zero, or 
proportional to the same terms, we also obtain the bounds on the second line of (3.11). 
Finally, since (see Section C) 

pAr,t||2, llPiV^tp, piV.tIb, \\'^xPN,t\\2, \\‘^xPN,t\\2, \\^xrN,t\\2 < 
we also obtain the estimates on the last line of (3.11). □ 


Next, we analyze the second term in (3.9). We have 

= I A) 

+ J dxdyN^^ViN^{x-y))^f{x)^fiy){s^,s^) 


+ Re 


c^) 


+ 


J dxdy N^^V{N^{x-y))ipf{x)(pf{y){si 
^ J dxdyN^^V{N>^{x-y)) [ipf^{x)ip^{y)alal + (pf {x)iff {y)a:,ay 


with given in (1.32). The properties of the generator are analyzed in the next 
proposition, which together with Proposition 3.3, implies Proposition 3.2. 

Proposition 3.5. Under the assumptions of Theorem 3.1, the operator d'^j^{t) defined 
in (1.32) satisfies the bounds 

< Cle^'''(A^+ 1), < Ce^'*'(A^+ 1)^ 

±[Ar,44W] ^ 1), ±[A^^4!'4^)] ^ 

44(^) ^ l4!3(0P < Ce^'''(A^ + 1)' 


24 


To prove Proposition 3.5, we will make use of bounds contained in the following 
lemma. 

Lemma 3.6. Let ji,j 2 S x M^). We consider the operators 

^ 1,1 = J dxdyN^^V{N^{x - y))ipi{x)ip2{y)a^{ji,x)a^{j2,y) 

^ 1,2 = J dxdyN‘^f^V{N^{x - y))ipi{x)ip2{y)a^{ji,x)ay 
^ 1,3 = J dxdyN^^V{N^{x - y))ipi{x)(p2{y)alay 


^ 2,1 = J dxdyN^f^V{N<^{x - y))ipi{x)ip2{x)a^{ji^y)a^{j2,y) 
^ 2,2 = J dxdyN^'^V{N^{x - y))ipi{x)ip2{x)a'^{ji^y)ay 
^ 2,3 = J dxdyN^'^V{N^{x - y))(pi{x)ip2{x)alay 


Then we have 

\{^|;,AJ,,^|;)\<C\\h\\2\\j2hM^ 

\{^P,Aj,2^P)\<C\\h\\2\MH4T2M^ (3.17) 

\{AAj,s^P)\ < C\\MH4T2Mm + 

for j = 1,2 and all ip & T. Moreover, 

— C'||jl||2lb’2||2ll7^l||^2||¥52|lH2 (-^ + 1)^ 

^j,2^j,2 + ^j,2^i,2 < C'||ji||2ll7’l||^2||<y92||^2 (A7 +1)^ (3.18) 

^j,3^p3 + ^j,3^i,3 < 11^2 1197211^2 (A7 + 1)^ 

for j = 1,2 and all ip £ T. 


Proof. We start with the bounds (3.17), for the case j = 1. We have 
\{ip,Ai^iip)\ < j dxdyN^^V{N^{x - y))\ipi{x)\\ip2iy)\\\a^{ji,x)fj\\\\a^{j2,y)'tp\\ 
< ||<y2i||oo||<y?2||oo j dxdyN^^V{N^{x - y))\\ji^c,\\2\\j2,y\\2\\{Af + 



<C||ji||2||j2||2||¥^l||H2||^2||H2||(A7+l)VVf 
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Moreover 


|(V',^1,2V')I < j dxdyN^^V{N^{x - y))\^pl{x)\\(p2{y)\\\a\jl^x)^l!\\\\ay 


< C'||(^i||oo||¥? 2| 


j dxdyN^dv{N^{x - y))\\ji,cc\\lm + 

dxdyN^^V {N^ {x — y))\\ay'il 


1/2 


1/2 




Finally, 


|(V’,^l,3'iA)| < J - y))|/3l(x)||v92(y)|||ayV'||||Oa:V’ll 

< ||</?i||oo||</?2||oo J - y))||ayV'|P 

The terms A 2 J, with j = 1,2,3 can be bounded similarly; we skip here the straightfor¬ 
ward details. We switch now to the estimates (3.18). Again, we consider the case j = 1. 
We have 

< J dxdydzdwN^^V{N^{x — y))N^^V{Nd[z — w)) 

X |/?l(x)||(/)2(y)||/9l(z)||(/>2(«;)|||o>*(jh,a:)o‘‘(j2,y)V’lll|o‘‘(jl,^)o‘‘(j2,i 

< C||(/?i||^||(/? 2 ||L j dxdydzdwN^^V{Nd{x-y))N^^V{NP{z-w)) 

X ||jl,x||2||j2,i/||2||jl,z||2||j2,«)||2||(AA-I- l)'i/’|P 
^ C'||(/?i 11^2 ||(/>2 11^2 II (AA -I- 1 ) 1 /)Ip 

X J dxdydzdwN^^V{N^{x - y))N^dv{Nd{z - wmh,,g\\j 2 ,, ' 

< C||jlf ||j2f||(/5l 111^2 ||(^2 11^2 ||(W+ l)V'f 

The contribution A^ lAi^i can be estimated exactly in the same way. Let us now consider 
the term Ai^sA^ 3. W^e find 

( 1 /), ApaA^ 3 ?/)) 

= J dxdydzdwN^^V{Nd{x - y))N^^V{Nd{z - w)) 

X ipi ( x ) ip2 { y )^ i { z ) if2 ( w ) ( V ’, alayala ^' ijj ) 

= J dxdydzdwN^^V{Nd{x - y))N^^V{Nd{z - w)) 

X (piix)ip2{y)^i{z)^2(w) (V’, alalaya^'i/j) 

+ J dxdydzN^^V{Nf^{x - y))N^f^V{Nd{z - y))ipi{x)\ip2{y)f'ipi{z){'ip,ala^'ip) 


|2 

!,u; II2 
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which leads us, by Cauchy-Schwarz, to ('i/’,^1,3^1 3V’) < C'||‘^i||^2||y^2||^2 
Since ^ 4 ^ 3 equals ^1,3, with ipi and ip2 exchanged, the operator 3^1,3 can be bounded 
in the same way. Also the estimates for the terms Ai^2A| 2, A| 2^1,2 and for the squares 
of the operators A2J, with j = 1, 2, 3 , can be obtained similarly. □ 

Proof of Proposition 3 . 5 . Writing cosh^j^ ^ = l+p7v,t) and recalling from Lemma C.l that 
sinhfcjy xM^), with a norm bounded uniformly in N, we notice that 

a sum of terms, each of them having the form of one of the operators Aj^i, 2; ^*,35 

for an i = 1,2. Since the solution of ( 2 . 10 ) satisfies ||(^^||j:/2 < the bounds 

( 3 . 17 ) in Lemma 3.6 immediately imply that 

and (since by the assumptions on ip we have ||(/9^||j:^2 < Ce^l^l) also that 

±[W, 45 W] < ± 4!2 < + 1). 

From ( 3 . 18 ), on the other hand, we find 

{62}? + if 

02,Nf +lf 

±[M\ 6 ^}{t)]<Ce^\^\{M+lf 


□ 


3.B Analysis of Tff ^g^f\t)TN,t 

We have 


TfG^Tt = 


1 

y/N 


j dxdyN^^V{Nf^{x-y))ipf{y) 


+ 


1 


^fN . 


X (a*(c)y) + a(4^))(a*(<) + a(<))(a(c)y) + a*(s)y)) 

dxdyN^^V{N^{x - y))(pf {y) 


X {o*{c^) + a{sf)){a{cf) + a*{sf)){a{cf) + a*{sf)) 
Writing the terms in normal order and decomposing 

{Sx,cf) = SN,t{y,x) + {sf,pf) = kN,t{y,x) + rN,t{y,x) + {sf,pf) 
we arrive at 


TfG^Tt = 


1 

Vn 


J dxdyN^V{N{x - 


y))Pt {y)kN,t{x,y)Tfa^Tt + I 1 .C. + 


( 3 . 19 ) 
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with 


^3,iv(i) = 


a/]V 



+ 


+ 


+ 


1 

y/N 

1 

y/N 

1 

y/N 


X - y))(Pt iy) 

X [«* (Cx )a* (Cy )a* (sf) + a* (c^ )a* (c^)a(cf) + a* (c^)a* (sf )a(s^ ) 
+ a* (Cy )a* (sf )a(sf) + a* (cf )a(s^)a(cf) + a* (c^)a(sf )a(c^) 
+a*(sf )a(sf )a(s^) + a{s^ )a{s^ )a(cf)] 

J dxdy N^I^V{N^{x - y))ipf{y){s^,s^){a*{c^) + a{s^)) 

J dxdyN^^V{N^{x-y))(p^{y) ) + fjv,t(x, y)] 

X (a(cf) + o*(sf)) 


J dxdyN^^V{N^{x-y))ip^ 


+ h.c. 


(3.20) 


The properties of are established in the next proposition. 

Proposition 3.7. Under the assumptions of Theorem 3.1, we find, for every <5 > 0, a 
constant > 0 such that 

±£3,N{t) < -^Vjv + CM^/N + Cse^\^\{M + 1) 

±[AA,T3,iv(i)] <6VN + CMyN + Cse^\^\{M+l) (3.21) 

±4,iv(t) < <5VAr + CM^/N + Qe^l*l(AA + 1) 

Moreover, we have 

|(V’i,^3,Af(i)V’2)| < ^{i-y)/2 [^^1’ {IC+N +l)fii) + {^2, + (A^T 1)^)V'2)] (3-22) 

for all ^/^i, '02 £ A". 


After applying the Cauchy-Schwarz inequality, we will estimate the cubic terms in 
£ 3 ,N{t) by quadratic and quartic terms. The quadratic terms can be controlled with 
Lemma 3.6. To bound the quartic terms, we will need the following three lemmas. 

Lemma 3.8. Let V € L^(M^), V >D and 0 < /3 < 1. Let ji,j 2 S L^(M^ x M^) with 


Mi := max 


sup 

X 


dy\ji{x,y)\^ 

y 


J dx\ji{x,y)\‘^ 


< oo 


for i = 1,2. Then we have 


J dxdyN^^V{N{x - y))\\a^{ji,x)aHh,i 


I dxdyN^^V{N{x - y))||a«(ji,,)a^0f < CMi\\{M + 1)0|| 


< Cmin(Mi||j2||2,Af2||ji||i)||(AA + 1)011^ 

2 
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for all ^ J-. These inequalities remain true if both operators applied to on the l.h.s. 
act on the same variable, i.e. 


J dxdyN^^V{N{x - y))\\a^{ji,x)a^{h,x)fj\\'^ < Cmin(Mi||j 2 ||i, M 2 ||ji|| 2 )||(AA + 1)V'||^ 
I dxdyN^l^V{N{x-y))\\aHji,x)axfjf <CMi\\{Af + l)fjf 


for all if ^ T. 

Proof. We observe that 


J dxdyN^^V{N{x - 


< 


aHji,x)aHj2,y)'ipf 
j dxdyN^^V{N{x - y))||ji,x||il|j2,y||ill(Ar + 1)V’II' 


< C'min(Mi||j2||i,M2||ji||^)||(AA+ 1)V'||^ 


and that 





The last two bounds are obtained similarly. □ 

In the next lemma, we control quartic terms where the arguments of creation and 
annihilation operators is the kernel cosh^j^^. 

Lemma 3.9. Let Q < jd < 1, V ^ > 0, and Vn be defined as in (1.36). Let 

kN,t be defined as in (2.11) and, as usual, let (y) = coshfc^j(y, x). Then, we have 

^ J dxdyN^^ViN^{x-y))a*{c^)a*{c^)a{c^)a{c(^) = VN + £{t) 

where, for every 6 > 0, we can find a constant Cs > 0 such that 

£{t) < SVn + Cs{M + if/N 

Proof. The lemma follows expanding a{c^) = Ox + a{Px) a{cy) = Oy + a{py) and 
applying Lemma 3.8 and Lemma C.l. □ 

Finally, we need a bound for the expectation of Vn in terms of the square of the 
kinetic energy operator. This is the content of the next lemma. 
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Lemma 3.10. Let V E > 0 and V be radially symmetric, 0 < /3 < 1. Then 


acr^ 1/0-7J 


x'-^x v y^y 


j dxdyN^^V{N^{x-y))\\axayilj\\‘^ < C j dxdy\\V^ 
Proof. We define 

Then we have Ag = V and therefore 

J dxdyN^^V{N^{x - y))\\a^ay'il;2f 


^ + C 


J dxdy\\V, 


x^x^y% 


(3.23) 


•V,, 


g{N^{x-y)) {axay'ip 2 ,axayif 2 ) 


= / dxdyN^g{N^{x - y)) [{Vyay'4^2-,axay'4^2) + {'^xO-xCiy, a^S/yay'if 2 ) + h.c.] 


< 


1/2 r 


1/2 


□ 


dxdy\\VxaxVyay'if2\\^ J dxdylgiNl^{x - y))\‘^\\axayfj2f 

+ iV^y dxdy\g{N^{x - y))\\\Vxaxayf; 2 f 
< C y dxdy\\VxttxVyayip 2 \\‘^ + C j dxdy\\Vxaxay'ijj 2 \\‘^ 

Here we used the fact that (from Newton’s theorem) \g{x)\ < C\x\~^. 

Proof of Proposition 3. 7. Let us decompose the operator {t) as 

= L{t) + Ci{t) + C2{t) 

with 

L{t) y dxdyN^^V{NP{x-y))y}^{y){s^,s^){a*{c^) + a{s^)) 

+ h.c. 

Ci{t) = ^ y dxdyN^^V{Nf^{x - y))^f (y) 

X [a*{Cx)o*{Sx)a{s^)+a*{c^)a*{s^)a{s^)+a*{cf)a{s^)a{c^) 
+a*{cy)a{s^)a{c^) + a*{s^)a{s^)a{s^) + a{s^)a{s^)a{c^)\ + h.c. 

C2{t) = ^J dxdyN^^V{NP{x-y))ipf{y) 


X (a* (cf )a* (c^)a* (sf) + a* (cf )a* (c^)o(c^)) + h.c. 


(3.24) 
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Let us begin with the linear terms in L(t). Using Cauchy-Schwarz inequality and 
Lemma 3.6 (together with the estimates in Appendix C for the Hilbert-Schmidt operators 
sinhfc^it is easy to check that 


{'ipi,L{t)'ip2) <CN ((V'l, (Ar+l)V'i) + (Ar +l)^/'2)) 


and, taking V’l = V’ 2 , that 

±L{t) < + 1) 

±[Af, L{t)] < ^ < CiV-V2e-^l*l(AA + 1) 

Next, we control the terms in C'i(t). Decomposing cosh^^^ = 1 + PN,t and applying 
Cauchy-Schwarz, it is clear that all contributions to {•ijji,Ci{t)ip 2 ) can be bounded by a 
sum of terms of the form 


1 

Vn _ 


I dxdyN^^V{N^{x - y))\ip^ 

J dxdyN^^V{Nl^{x - y))\\a^{j 2 ,x)a{h,y)'ip 2 \ 


1/2 


Vn 


J dxdyN^^V{N^{x - y))\(p^{y)\‘^\\J{ji^^)'ilji\ 


-I 1/2 


1 

_ 


j dxdyN^^V{N^{x - y))||o**(j2,x)ayi/'2|P 
j dxdyN^^V{N^{x - y))\ip^{y)\‘^\\aoo'tpi\f 


1/2 


Vn 


J dxdyN^^V{N^{x-y))\ipf^{y)\^\\aM\^ 


dxdyN^^V{N^{x - y))||a**(j2,x)a**(j3,?/)'02"^ 
-11/2 


-|l/2 


X 


dxdyN^^V{N'^{x - y))\\a^ {j 2 ,x)ayilJ 2 f 


1/2 


with appropriate ji, j 2 , js £ L^(M^ x M^) satisfying ||ii ||2 < and 


Mi = max 


sup / dy|ji(x,y)|Usup / dx\ji{x,y)\^ 

. X J y J 




for all i = 1,2, 3. Using Lemma 3.8, we conclude therefore that 

|(iAi,C'i(t)V^2)| < CiV-V2e^l*l||(A/ + l)V2^,||||(A/+l)^2|| 

< [(V'l, {Af + l)V'i) + {ip2, {Af + 1)^V’2)] 


(3.25) 
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Taking -01 = '02 the first inequality shows that, for all <5 > 0, there exists a constant 
Cs > 0 with 


±Ci{t) < dM^/N + +1) 

±[AA,Ci(t)] <6AfyN + Cse^^^^{Af+l) 

±Ci{t) < + Cse^^^^{J\f +1) 

Finally, we study the term 6*2 (t). Using Cauchy-Schwarz, we obtain 

^ j dxdyN^/^V{N^{x-y))ip^{y){'ilj,a*{c^)a*{c^)a{c^)'ip) 
j dxdyN'^^V{Nl^{x-y))\\a{c^)a{c^)i)f 
+ ^l dxdyN^f^V{NPix-y))\^f{y)\^\\a{c^W 
for every k > 0. With Lemma 3.9, we find, for every 5 > 0, a constant > 0 such that 
/ dxdyN^dv{Nd{x -y)) [(^f (y)a*(c^)a*(c^)o(cf) + h.c.] 

< 5 Vjv + CAf^/N + {Af + 1) 

Similarly, one can show that, for all <5 > 0, there exists > 0 such that 

/ dxdyN'^dy(^j^0(^^ [p^{y)a*{c^)a*{c^)a*{s^)+ h..c] 

< ^Vat + CAf‘^/N + Qe^l*l(W + 1) 

We conclude that, for all <5 > 0 there exists C 5 > 0 with 

±C2{t) < 5Vn + CAf^/N + C5e^\^\{Af + 1) 

± [Af, C2{t)] < 5Vn + CAP/N + Cse^\^\Af + 1) 

±C2{t) < <5Viv + CAf‘^/N + C5e^\^\{Af + 1) 


To prove (3.22) for the term C 2 {t), we need to proceed differently. We observe that, 
for any , 02 S -U, 



„N\ 


(3.26) 


X - y))^f (y) (01, a* (c)) )a* (p^ )a(c)y )02) 
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To bound the second term, we observe that 


J dxdyN^^V{Nf^{x-y))ip^{y){'ilJi,a*{c^)a*{p^)a{c^) 2 p 2 ) 

< J dxdyN^^V{Nl^{x-y))\^pf^{y)\\\p^\\ 2 \\a{c^)^pl\\\\a{c^){^^ + l)^^'^^p 2 

< ^ + (t/' 2, (AA + l)^V’2)] 


As for the first term on the r.h.s. of (3.26), we recall the definition (3.23) of the function 
g with /\g = V. Putting v = Vgf, we find that V ■ v = V. Hence 


1 


^ J dxdyN^I^V{N^{x-y))ipf{y){'ipi,a*{c^)ala{c^)ijj 2 ) 

dxdyN'^^Vy \v{N^ix-y))] (y)(V’i,a*(c^)o*a(cf)V' 2 ) 


Vn . 
7^, 


+ 


1 

Vn 


dxdyN'^^v{N^{x - y)) ■ (?/)(?/;i, a*(c^)a*a(c^)V’ 2 ) 

[ dxdyiff {y)N'^^v{Nf^{x - y)) ■ {a{c^)Vyay'ijji,a{c^)'ijj 2 ) 


Therefore 
1 


^ J dxdyN^^V{Nf^{x - y))ipf{y){'ipi,a*{c^ )a*a(cf )'02) 

/ dxdyN‘^^\v{N^{x-y))\\\a{c^){Af + iy^‘^'ilj2\\ 

|V(/?f(y)|||a(c^)(AA + l)-^/2ayV^i|| + \ipf^{y)\\\a{c^)iM + ly^^Wyayil^iW 

(3.27) 


The contribution of the first term in the parenthesis can be estimated by 


j dxdyN‘^^\v{Nf^{x - y))\\Vipf {y)\\\a{c^){Af + 1) ||a(cf )(AA + l)^/2^2 

1/2 


J dxdy\\a{c^){Af +1) 

j dxdy\v{Nf^{x - ?/))|^||a(cf )(A/+ l)^/^'02 


1/2 




Here we used the bounds ||u||oo < C'll^lls ||i ;||2 < <1/111^116/5) Hi^'t c™ b® proven with 
the help of the Hardy-Littlewood-Sobolev inequality. 
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As for the second term in the parenthesis on the r.h.s. of (3.27), we have 

j dxdyN‘^^\v{Nf^ix-y))\\ipf{y)\\\a{c^){^f+l)~^^‘^Vyay^pl\\\\aic^){^f+l)~'^/'^'^|J2\\ 

r f -11/2 


J dxdy\\a{c^){Af + 1) II" 

X J dxdy\v{N^{x - y))\‘^\\a{c^){J\f + 
< Ce^l*liV^/2 [(V-i, /C^i) + {^2, + l?^2)] 


1/2 


We find that 
1 


^ J dxdyN^^V{N^{x - y))cpf {y) (i/^i, a* (cf )a* (c^)a(cf) 1 /^ 2 ) 

< C7e^|tlAr(/3-i)/2 (/C + AA + I)!/-!) + (i/>2, (AA + 1)2V;2)] 


(3.28) 


Since the other term entering in the definition of C 2 {t) in (3.24) can be bounded similarly 
and the hermitian conjugates can be treated using Lemma 3.10, we conclude that 


KV’I, C2(t)V>2)| < {1C+M+ i)V’i) + ( 1 A 2 , {1C+M + if ^ 2 )] 


□ 


3.C Analysis of 

Finally, we consider the conjugation of the quartic terms. Expanding the products, and 
writing all terms in normal order, we obtain 

2T:Qf{t)Tt 

= ^j dxdyN^^V{N^{x-y)) [\{sf + \{s^, sf)\'^ + {sf, sf ){s^, sf)] 

+ ^ / dxdyN^^V{N/^{x-y)) [{cf ,sf)a*{cf)a*{cf) + {s^,cf)a{cf)a{cf)] 


+ 2V7V + £A,Nit) 


(3.29) 


where we defined with 

= J dxdyN‘^V{N{x - y)) 


2{sf ,sf)a*{c^)a*{sf) + 2{sf ,sf)a{s^^)a{c^^) + 2{s^^ ,sf)a*{cf)a*{s^^) 


„N N 


N\ 


„N\ 


W „N\*/N\ 


^N\ 


+ 2(sf,s^)a(s(;)a(ciy) + 2(s(;,sna*(4^)a(c^^) + 2(s(;,s^^)a*(c^^)a(c(;) 


r,N\ 


^N\ 


r.N „N\*r„N 


.N\ 


r,N „N\*f„N\ 


^N\ 


+ 2(s-,s-)a*(s-)a(s-) + 2(s;\s-)a*(s-)a(s-) + (c;\s-)a*(c-)a(s-) 
+ {sf,cf)a*{sf)a{c^) + {sf ,cf)a*{sf)a*{s^) + {cf ,s^)a{sf)a{sf) 


y 


+ {sf,cf)afsf)a{cf) + {cf, sf)a*{cf)a{s^) 
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and 


(3.30) 




Here, the term 


[ dxdyN^f^ViN^ix-y)) 


X a* {Px )a* [cy )a{c^)a{c^) + a^a* {p^)a{c^ )a(c^) 

+ alala{p^ )a{c^) + ala*yaya{p^] 

contains the contributions arising from a*{c^)a*{cy)a{cy)a{c^), after removing Vat, 
while 

+ a* (c^ )a* [c^ )a* (sf )a{c^) + a* (cf )a* [s^ )a* {s^ )a{s^) + a* (cf )a* [s^ )a{s^ )a(cf) 
+ a* (c^ )a* {s^ )a{s^)a{c^) + a* (c^ )a{s^ )a{c^)a{c^) + a* (c^ )a* (s^ )a* (s^ )a{s^) 

+ a*{Cy )a* {s^ )a(sf )a(cf) + a* (c^ )a* {s^)a{s^)a{c^) + a* (c^ )a{s^)a{c^ )a(cf) 

+ a*{s^ )a*{s^ )a{s^ )a{s^ ) + a*{s^ )a{s^ )a{s^ )a{c^ ) + a*{s^ )a{s^ )a{s^ )a{c^) 


+ a{sx)a{sy)a{c^)a{c^) 


(3.31) 


Proposition 3.11. Under the assumptions of Theorem 3.1, there exists, for all 5 > 0, 
a constant Cs > 0 such that 

±£4,N{t) < SVn + {Af + if/N + {Af + 1) 

± [Af,£4,Nit)] < 6 Vn + Cse^\^\iAf + if/N + CsN-^+^e^^^^iAf + 1) 

±£4{t) < SVn + {Af + if/N + {Af + 1) 


Furthermore, 


\{ifi,£4,Nm2)\ < {{Af + lf/N)f,f 


+ {ipi,{Af + ffji) + {'if2,{IC+Af + if 1 ^ 2 ) 


(3.32) 


for all V'l, i /’2 £ !£■ 

Proof. We start with the terms in f (t). We write 

£^^l{t) = £i%\t) + Si^fft) 
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with 


(t) J dxdy N^^V{N^{x - y)) 


2{s^,s^)a*{c^)a*{s^) + 2{s^,s^)a{s^)a{c^) + 2{s^,s^)a*{c^)a*{s^) 


+ 2(s^, s^)a{s!^ )a{c‘^) + 2{s !^, s" )a*(c" )o(4'') + 2(s", s‘^)a*{c'^ )a{c'^) 




„N\ 


y 

N\*fN 


.N\ 


„N N\*t„N\ 


„N\ 


+ 2(s^,s^)a*(s^)a(4'^) + 2(4\4'')a*(4'')a(4'') 




-'y ) -y 

„N\*r„N\ 




I dxdy N^f^V{N^{x - y)) {c^, s^)a*{c^)a{s^) 


N 


+ {sx ,Cy) a*(s^vMcx) + (4\<'')«*«)a*(4'') + 4'')o(sf )a«) 






-'y 


V 


+ (44<)a4444<) + «,44^44 )«(44 

Using Lemma 3.6 and the bounds for the Hilbert-Schmidt operator sinh^j^j, we easily 
find that 

(41,45/^(442) < C'a^"4^i*I||(aa +i)^/^i/;i||||(AA +i)^/2V'2|| 

< [(4i, {J\f + i)4i) + (42, {J\f + 1)42)] 

Taking 'ipi = 42, we immediately obtain 

±45/54 < ca^"4'^i*I(aa + 1) 

± [aa, 44^541 < C’a^“4^i*I(aa + 1) 

±4^’^54 < ca^"4^i*I(aa + 1) 

(2 2 ) 

Let us now consider the terms in They all contain the inner product 

{c^,s^) = kN,t{x,y)+rN,t{x,y) + {p^,s^). Since \rN,t{y^^)\ and |(p^,s^)| are bounded, 
uniformly in N, the contributions arising from these terms can be dealt with, as we did 
above, for the term T 4 (t). For the contributions proportional to kj\f^t{x, y), we use the 
bound \kN^t{x,y)\ < C\x — y\~^. Paying the price of an additional factor we can 
deal with these terms as we did above, just replacing the potential V{x) with V{x)/\x\ 
(which is still integrable, by assumption). We conclude that 

I( 41 , 445442 )I < + i)^/2V’illll(AA + i)^/^42|| 

< [(4i, (A^± i)4i) ± (42, (A^± 1)42)] 

As usual, we also obtain 

±4454 < C'Af-(^-^)e^l‘l(W+ 1) 

W,4454] <<TAi-(^-^)e^l‘l(W+l) 

±4^’^54 < + 1) 
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Finally, we consider the quartic terms in Recall the decomposition (3.30). 

From Lemma 3.9, we find for every (5 > 0 a constant (7^ > 0 such that 

±£i%\t) < 6 Vn + C5e^l*l(Ar+ l)ViV 

< -5Viv + Cse^\^\iAf+lf/N 

(notice that [AA,= 0). In order to show that satishes the estimate (3.32), we 
remark that, from Lemma 3.10, 

\{'llJl,VN'lp 2 )\ < (V'l, 

<VatV'i) + ;^ J dxdyN^^V{Nl^{x-y)){il^ 2 ,ala*yaya^'il; 2 ) 

< [(V'l, VatV'i) + 0^2, (/c +AA+ 1)^2)] 

Analogously, with Lemma 3.9, we find 





2 ^ 

{lp2,yNlp2) H-^^ll(AA+ 1)V'2|P 


< ^ VjvV’i) + j^WiM + + (V'2 ,(a: +AA + i)2 V’2) 

The last two equations show that 

|(V’i,^l^jV^^(t)V’2)l < (V’l, VatV’i) + ^||(AA + l)V^if + (V^2, (/C+AA + 1)^2) 

We switch now to defined in (3.31), which can be further decomposed as 


^4,Af^ (^) ~ + ^2(i) + ^3{t) 


with 


Ai(t) = ^ / -y))«*(s^)«*('S^)a(sf)a(s^) 

^ J dxdyN^^V{N^{x-y)) 


+ 


X 


a* (cf )a* {s^ )a* (s^ )a(s") + a* (c" )a* (s" )a* (s" )a(s^'') 


* / ^N\ 




+ a* )«( 4 '' )a{sy )a(c)y) + a* (s" )a(s)y )o(s 7 ' )o(c 7 ') 


„N\ 


„N\ 


„N 


*/„N\ 


„N\ 


„N\ 




+ 


X 


^ J dxdyN^^V{N^{x-y)) 


a* (cf )a* {s^ )a(s" )a(c^'') + a* )a* (s" )a(s()' )a{£^) 


y 


„N 


* („N\„*f „N\ 


y 


-'y 


+ a* {Cy )a* {s^)a{s^)a{c^) +a* {c^)a* (sf )o(sf )a{c^) 
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and 


^2(t) = ^ J dxdyN^^V{Nf^{x-y)) 

X a* {cx )a* (cy )a* {s^)a* (sf) + a(s^ )a{s^)a{c^ )o(cf) 

^ 3 (i) = ^ J dxdyN^^V{N^{x-y)) 

X \a* (c^)a* (cf )a* (sf )a(c^) + a* (cf )a* (c^)o* (sf )a(c^) 


+ a*(c^ )a« )a(c''' )a(c") + a* {c!^ )a(4'' )a(c" )a(4) 






„N 


*t„N\ 


„N\ 


„N\ 


In Ai{t), we collected terms with 0, 1, and some of those with 2 kernels cosh^^ ^ (namely 
the terms with 2 kernels that can be separated by Cauchy-Schwarz). In A 2 {t) we collected 
all other terms with 2 factors cosh^j^^. In ^ 3 (t), we included all terms with 3 factors 
coshfcj^j. Terms in Ai{t) can be controlled with Lemma 3.8. We hnd 






-WiM + 1)M" + m + l)i;2f 


Taking ipi = '4^2-, we get 


±Ai{t) < + if/N 

±[AA,Ai(t)] < Ce^l*l(AA+l)ViV 
±ii(t) < Ce^l*l(AA+ l)ViV 

Next, we consider terms in A 2 {t). By Cauchy-Schwarz and Lemma 3.9, we have 

^ J dxdyN^f^V{Nl^{x-y)){4}i,a*{cf)a*{cf)a*{sf)a*{sf)4)2) 



(3.33) 


5N 




Hence, choosing 6 = N we hnd 
1 


N 


/ dxdyN^^V{N^{x - y)) (V'l, a*{cf )a* [cf )a* {sf )a* {sf )V’ 2 ) 


C7e^ld 


< 


y/N 


{4ji,Vn'iPi) + —||(AA + l)4f -k (V’2,(AA + 1)^2) 
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The hermitian conjugated term can be bounded using Lemma 3.10; we get 

^ J dxdyN^<^V{Nf^{x - y)){ipua{s^)a{s^)a{c^)a{c^)Tp2) 


Hence, we find 

{{-Ipl, A2{t)'lp2)\ < 


r 1 

^ (tAi, VatV’i) + — ||(AA + l)V’if + (t/'2, (/C + AA+ 1)^2) 


Going back to (3.33) and choosing V’l = V’ 2 ) we obtain that for every 5 > 0, there exists 
6*5 > 0 with 

±A 2 {t) < 6 Vn + G 5 e^l*l(AA + if/N 
±[M,A2{t)] <6VN + Cse^^^\{M+lf/N 
±A 2 {t) < 6 Vn + G 5 e^l*l(AA + if/N 

Finally, we bound the term A 3 . To this end, we observe that 


- j dxdyN^^V{N^{x-y)){'ipi,a*{c^)a*{c^)a*{s^)a{c^)'ilj2) 

dxdyN^^V{N^{x-y))\\a{c^)a{c^)'ilJif 

<^/ dxdyN^i^V{Nl^{x - y))\\axay'tpi\\‘^ 

+ + ^^\m+1)4,2 f 

Similarly, the complex conjugated term can be controlled by 

^ J dxdyN^^V{Nf^{x-y)){4}i,a*{c^)a{s^)a{c^)a{c^)4}2) 

npK\A 

< 5 {4:2, Vn4^2) H- 77 —\\{Af + 1 )V’ 2 |P H-7777—II(AT + l)^/’i|P 


N 


SN 


Hence, we obtain that, for every (5i, <52 > 0, 

(V’l, A 3 (t)V’ 2 ) -^y dxdyN^^V{N^{x - y))\\axay 4 ;i\f 


+ 


N 
Ce^l^l 


+ ^y dxdyN^^V{N^{x - y))\\axay 4 ; 2 \\‘' 


N 


1 \ 

62 + y] ||(AA+l)^ 2 f + 


5i + 


52 


ll(AA+l)V’i|P 
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Choosing 6i = N and 62 = we find 


(^l,^3W'02) 


< 


C7e^|i| 

y/N 


(V'l, VtvV'i) + —{'ipi,{M+ l)Vi) 


(V'2,(/C+AA+1)V2) 


If instead we set di = <52 = <5 < 1, we obtain, with ipi = 1 ^ 2 , that for every <5 > 0 there 
exists > 0 with 


±A3{t) < 6 Vn + Qe^l*l(Ar+ l)ViV 
± [M,A3{t)] <5 Vn + Cse^\^\{M+l)^/N 
±A3{t) <6 Vn + Cse^^^\{M+lf/N 


□ 


3.D Analysis of {idtT^^)TN,t 
We set 

B{t) = ^ J dxdy(^kN,t{x,y)ala*y - kN,t{x,y)a^ay 

and 

-B(t) = ^ J dxdy[kN,t{x,y)ala*y-kN,t{x.,y)axay^ 

Then Tfq^t = exp(S). Defining ad^(C') = C and ad^’^^(C') = [i?,adg(C')] we have 

/■i _ 

{idtT*^,t)TN,t = - I = J] i-^adS(i?) (3.34) 


n>0 


on the domain T>{M) of the number of particle operator (this can be shown as in [6, 
Lemma 2.3]). In the next lemma, we collect some bounds for the operators ad^(i?), 
whose proof can be found in [6, Lemma 6.9]. 

Lemma 3.12. For each n G N there exist /n,ij/n ,2 £ L^(M^ x M^) such that 


j 5 


where 


a(l^{B) = ]^ J dxdy(^fn,i{x,y)alal + fn,2{x,y) 
a(f^{B) = ^ J dxdy(^fn,i{x,y)a*ay + fn,2{x,y)axal 

ll/n,i||L2xL2 < 2”llA:Ar,tll2 ||^Af,t||2 


for all n even and 
for all odd n , (3.35) 

(3.36) 


for all n >0 and i = 1,2, 


u 


n,i\\2 ^ 


\\kN,th ifn = 0 

4 "llA:Ar,f II2 ^ (||^Af,t||2||^Af,tII2 + ||^Af,t|l2) */ra > 1 , 


(3.37) 
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and 


j dx\fn,i{x,x)\ < 2 ”||A;Ar,i|| 2 ||feAr,t ||2 , 

j dx\fn,iix,x)\ < 4”||A;Ar,i||2"^ (||%,till + IlfcAT.tlbll^Ar,til) 

for all n > 1. 


(3.38) 


Using the bounds from the last lemma, we can now control the operator [idtTf^ 
appearing in the generator C^it). 


Proposition 3.13. Under the assumptions of Theorem 3.1, we have 

±{idtTf,^t)TN,t < + 1), + 1)' 

± [Af, {idtT;,^t)TN,t] < + 1), ±[A^^ {idtT*N,t)TN,t] < + 1)' 

±dt[{idtTf,^t)TN,t] < Ce^l'l(AA + 1), [dt[{idtTf,^t)TN,t]Y < Ce^^^^{Af + 1)' 

Proof. For any /i ,/2 G x M^) with f dx|/ 2 (x,x)| < oo we have 


{ip, J dxdy{fi{x,y)a*al + f 2 ix,y)ao,ay)ip) < (H/ilb + ||/ 2 || 2 )(V’, (A^+1)V') 


and 


{f;, J dxdy{fi{x,y)alay + f 2 {x,y)a,r^a*y)'il)) 

< (II/ 1 II 2 + ||/ 2 || 2 )(V',Art^) + j dx\f 2 i. 

From Lemma 3.12 and (3.34), we find therefore 


n>0 ^ 

n>0 ^ n>l 

<e2|l"^.‘IH|fcw,t||2||(AA+l)V2^||2 
<Ce^l*l||(AA + l)V2^||2 


(3.39) 


(2||%,t||2)^ 

(2n)! 


llWlbll'^ll 


(3.40) 


The bound 

\{f;,[Af,{idtT*^^,)TMM)\ < Ce^'*l||(AA + l)V2^f 

follows similarly, since the commutator with Af is just eliminating the contributions 
ad^(i?) for all odd n. 
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To prove the bound for , we use again (3.34): 

EE („(n| fa'ni (V',ad^(.B)ad^(.B)V^) 

n>0 ^ ^ '' 

SE 


n>0 m>0 

|| adS(i?)^|| ||adg(i^)V;|| 
(n + 1)! (m + 1)! 


(3.41) 


We claim that 

||ad^(S)V>f = (V’, ad^(i?) >) <C(||/„,i||2 + ||/n,2||2)'(V’,(AA + l)2V>) (3.42) 


for n even and 


||ad^(i?)V’r = U, adS(i?) 


■ l2 


V’ 


< (||/n,l||2 + ||/n,2||2)' (V’, (AA + 1)^) + C 


for n odd. The bound for n even follows because 


J dx\fn,2{x,x)\ 


(3.43) 


{'Ip, ad%{B) -Ip) = j dxdydx'dy' fn,i{x,y)fn,i{xPy'){'ip,alalal,a*y,'ip) 

+ J dxdydx'dy fna{x,y)fn, 2 {x',y') {ip, a^aya^iay'Ip) 


(3.44) 


+ J dxdydx'dy' fn,iix,y)fn, 2 {x',y')ala*ya^'ay' 

+ j dxdydx'dy' fn, 2 {x,y)fn,iix',y'){ip,axayal,a*y,'ip) 

The absolute value of the first term is bounded by 
J dxdx' Wa^a^'ipWWa* {fn,i,x)a* {fn,i,x')'^\\ 

<C\\{J\f+l)'lp\\ J dxdx'\\fn,l,x\\ 2 \\fn,l,x'\\ 2 \\axax''tp\\ 

The contribution of the second term on the r.h.s. of (3.44) can be estimated similarly. 
The third term, on the other hand, is bounded by 


J dxdydx'dy'\fn,iix, y)\\fn, 2 ix', y') 


< J dxdydx'dy' [\fn,i{x,y)\'^\\ax'ay''ip\\‘^ + \fn, 2 {x',y')\'^\\axay 
<{\\fn,lf2 + \\fn,2f)\m+l)^Pf 
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To control the last term on the r.h.s. of (3.44), we observe that 
J dxdydx'dy’fn^ix, y)fn,i{x', y'){il}, a^ayal,a*y,'il^) 

= J dxdydx'dy'fn, 2 {x,y)fn,iix',y'){a^'i’,a::cayap'ijj) 

+ J dxdydy'fn, 2 {x,y)fn,iix,y'){iJ,ayap'ilj) 

+ J dxdydy'fn, 2 {x,y)fn,iiy,y'){ilJ,axap'ilj) 

= j dxdydx fn,2{x,y){a:c'il;,axaya*{fn,i,x')i>) 

+ J dxdydy'fn,2ix,y)fn,iix,y'){'4’,apay't{j) + 

+ J dxdydy'fn, 2 ix,y)fn,i{y,y'){'ilJ,apax'ip) + 

This implies that 

f dxdydx'dy'fn,2{x, y)fn,i{x', y'){ip, a^ayapapif;) 

= j dxdydx'lfn,2i^,y)\\\ax'{J^ + l)^^‘^'ilj\\\\axaya*{fn,i,x')iJ^ + 

+ j dxdydy'\fn, 2 {x,y)\\fn,li^,y')\\\ay'tp\\\\ay>'tl;\\ + ||l/’|P||/n,l|| 2 ||/n, 2||2 

+ j dxdydy'\fn,2{x,y)\\fn,liy^y')\\\0‘y''^\\hx'tp\\ + llV^f ll/n,2||2||/n,l||2 

< C'(ll/n,l||2 + ll/n,2||2)ll(-^ + 


J dxdyfn,iix,y)fn,2ix,y) 

J dxdyfn,2{x,y)fn,iiy^^) 


and concludes the proof of (3.42). Eq. (3.43) can be shown analogously. Combined, 
they imply that 

and thus, inserting in (3.41), that 

[{idtTp^t)TN,t]^ +if 

Let us now consider the bound for the commutator [N"^, {idtTp ^)T]\f^t] in Proposition 
3.13. We observe that 

(</>, [N \= !(</.,V[V,(!8,rj_,)r„,,]./.) + (</.,|V,(!8,rjj)rjv.,]Vv-)| 


n even 

scE 

n even 


1 


(n + 1)! 

1 


(n + 1)! 


(||/n,l||2 + ||/n,2||2)||(AA+l)-"2 
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Finally, we observe that the bounds involving the time-derivative of {idtT^ in 

Proposition 3.13 can be proven similarly, taking the time derivative of the expressions 
for adg(i?), and using the bounds for ||/n,i ||2 and f dx\fn,iix,x)\ in (3.37) and (3.38) 
and, finally, using the estimate for ||A:jv,t ||2 proven in Appendix C. □ 


3.E Proof of Theorem 3.1 

Finally, we combine all results of the previous subsections to prove Theorem 3.1. 
From (3.5), (3.6), Prop. 3.2, (3.19), (3.29) and (3.34), we conclude that 

^N{t) = 1lN{t) + {idtT^-i-)T]\[^t + + £4,N{t) 


with 


+ 


+ 


Ai = + h.c.] Tat,* 

+ j dxdyN^^V{N^{x -y)) {y)kN,tix,y)T^^i-alTN,t+ 

A 2 = / dxdy {-A^kN,tix,y))a*xal + / dxdy {-A^kN,tix,y))axay 

J J (3.46) 

^ J dxdyN'^^V{Nf^{x-y)) [ip^{x)^p^{y)alal + (p^{x)(p^{y)axay\ 

^ J dxdyN^^V{N^{x-y)) 

X [{Cy,s^)a*{c^)a*{c^) + {s^,c^)a{c^)a{c^)] 

We observe that, by the definition of kj^^t, 

Ai = Vn j dxdyN^dv[Nd{x - y))ipf {y)uN,£ix - y) 

X [^f (y)(^f (x) - (pt{{x + y)/2f] T^a^^Tt + h.c. 
Using the bounds in Lemma 2.1, we find 

I dyN^^V{Nd{. - y))p>^{y)uJNA- “ iv) [^f (y)^f (•) “ ((• + y)/2)^] ^ 

and therefore 

|(V’i,AiV’ 2)1 <CiV-V2gi^M||(A7 + l)V2^i||||(AA + l)V2^2|| (3.47) 
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Let us now consider the second term in (3.46). Using the permutation symmetry of 
T, we find 


/ 


dxdy{A^kN,t)ix, y)ala* 

= j dxdy{AykN,t)ix,y)ala* = ^ j dxdy[{Aa, + Ay)kN,t{x,y)] 


dxdy 


‘^Ax—y + ^A(^x+y)/2 ) kN,t{x,y) 


= -^^/ dxdy{AojN,d{x-y){(p^ {{x + y)/2)f ala. 

N f 

- y / dxdyujN,(.{x - y) 


* * 
axay 


* * 
Ojx^y 


(3.48) 


X [v>H{x + y)/2)ALp^{{x + y)/2) + ((x + y)/2)p] a*o* 


Moreover, we have 

^ J dxdyN^^V{Nd{x-y))Lp^{x)ip^{y)alal 


i I dxdyN^^V{Nf^ix-y))ipf^i{x + y)/ 2 fala*y + Bi 


(3.49) 


where 


Bi = -- J dxdyN^^V{Nf^{x-y)) [ip^{{x + y)/2f - ipf {x)ipf {y)] 


* * 
Qj^ dy 


and 


^ I dxdyN^^ViNf^ix - y)){c^,s^)a*{c^)a*ic^) 


with 


B 2 = 


= ~\j dxdyN^^V{Nd{x-y))ujNAx-y)iV’^^ii^ + y)/‘^)fo.*x(^*y + B2 

y J dxdyN^^V{N^{x - y))[{p^,s^) + rN,t{y,x)]a*{c^)a*{c^) 

dxdyN^'^V{N^{x-y))ujN,i{x-y){y^^{{x + y)/ 2 )fa*{j)^)a*y 

dxdyN^^V{N^{x-y))u!N,i{x-y){(pf{ix + y)/2))‘^a*{c^)a*{p^) 

By the definition of = 1 “ ^N,t, we conclude that 
N 


^2 = y y dxdyujM,i{x - y) 

X [^t{{x + y)/2)Aip^{{x + y)/2) + ((x + y)/2)\^] a*ay 

+ NXN,e J dxdyl{\x-y\ < £){ip^{{x + y)/2))" 

+ Bi + B 2 + B 3 + h.c. 


(3.50) 


\ 2 ^ ^ * 
dj. dy 
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with 


B 3 = -NXn^i j dxdyujN,iix - ((x + y)/ 2 )fala*y 

Next, we show that Bi,B 2 ,B 3 are small. To control Bi, we define f^v,* through 

k - y\ y) ■= ‘ft{x)^tiy) - + y)/‘^) 

We set W{x) = V{x)\x\, we define 


sw = -i / 


dy- 


x-y\ 


W{y) 


and V = 'Vg-, then we have V • u = W. Integrating by parts, we find 


\{'lpl,Bi'lp2)\ = 


2NP 


J dxdyN^^V{N>^{x - y))iV^|x - y\ ft,Ar(x, y)(V'i, V' 2 > 


< 


j dxllaa^-^ill ||a*(iV^^u(iV^(x-•))Va:ft,Ar(x,-))^2|| 

dxWVxaxipiW \\a*{N'^^v{N/^{x - ■))ft,N{x,-))ip2\ 


2NP 


(3.51) 


Since ||?;||2 < C'llMl^lle/s we obtain 

\\N‘^f^v{Nf^{x - y))ft,N{x,y )\\2 < 

\\N^^v{N^{x - y))Vxft,N{x,y)h < 

From (3.51), we conclude that 

KV’i, + Bl)i;2)\ < CiV-^/2gi^|t| (M + JC + l)V.i) + (^ 2 , (AA + /C + 1)V’2)] 


To bound B 2 , we notice first that 

^ J dxdyN^^V{NP{x-y)) [{p^, s^) + rN,t{y]x)] (V^i, a*(cf )a*(c^)V> 2 ) 

npK\t\ 

- jYl/2 [^iV + {M + if /N] V’l) + (V’2,V'2)] 

as it easily follows from Lemma 3.9. Moreover, we find 

i J dxdyN^f^V{N^{x-y))ujN,fx-y){ipf{{x + y)/2)f{i)i,a*{pf)a*yi)2) 

pK\t\ r 1 

- J (^^dyN^f^V{N^{x - y)) Wpfh ||a^V^i||||(W + l )^/^^2 

< \\{Af + if/^M II (AA + 1 )^/^V '2 II 
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and an analogous estimate for the term proportional to a*{c^)a*{py). Hence, we con¬ 
clude that 


< Cmax ^K\t\ 


X [(V-i, (V^ + (AA + 1)VA^ + l)V^i) + (V^2, (Viv + {^f+ l?/N + 1 )^ 2 )] 


Finally, we control H 3 . From the bounds in Lemma 2.1, we immediately find that 


KV'l, (Hs + H3*)V^2)| < CN-^ [(V'l, {J\f+l)'iPl) + (^2, (Ar+ 1)V’2)] 
Combining (3.45) with (3.47), (3.50) and with the bounds for Bi,B 2 ,B^, we find 
C,N{t) = r]iq{t) + C2,N{t) + Vat + £N{t) 


with 


C2,N{t) = {idtT;r,t)TN,t + C?i,{t) 

N f 

-h y / dxdyixiN,(.{x - y) 

X + y)/2)Apf^{{x + y)/2) + \Vip^{{x + y)/2)p)a*a* -h h.c.] 



(3.52) 


and where £N{t) = £N, 3 {t) + £N, 4 {t) + Ai + Bi + B 2 + B^ satisfies the estimates (3.2), 
(3.3) in Theorem 3.1. 

To complete the proof of Theorem 3.1, we still have to show that the last two terms 
on the r.h.s. of (3.52) satisfy the bounds (3.4). But this fact follows easily, since, from 
Lemma 2 . 1 , Nujn/{x — y) < Cldx — y\ < £)/\x — y\ is square-integrable and < C 

is of order one. 

4 Growth of Number of Particles and Energy 

We apply the estimates of the Section 3 to show a bound for the growth of the expectation 
of the number of particles with respect to the fluctuation dynamics generated by CN{t). 
Moreover, we prove bounds for the growth of the expectation of the number of particles, 
of the kinetic energy and of their square with respect to the dynamics generated by the 
quadratic part of the generator £2 

We need, first of all, to compare kinetic and potential energy operators with the 
generators CN{t) of the full fluctuation dynamics and with its quadratic part £ 2 ,iV- The 
following proposition follows directly from Theorem 3.1. 
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Proposition 4.1. Under the same assumptions as in Proposition 2.2, there exists a 
constant C > 0 such that 

i(/C + Viv)-Ce^l'l [Af + l+Af^N] <CN{t)-VN{t) . 

2. (4.1) 

< 2(/C + Vtv) + [AA + 1 + AA^/iV] 

Moreover, we have 

/C-Ce^l‘l(AA+l) <C2,N{t) </C + (Ae^l*l(AA + l) (4.2) 

and 

-C'e'^l^l(Ar+l)2 <Cl^{t) <2/C2 + Ce^l*l(Ar+l)2 (4.3) 

Next, we control the growth of the expectation of the number of particles and of the 
energy with respect to the fluctuation dynamics UN{t\ s). 

Theorem 4.2. Under the same assumptions as in Proposition 2.2, there exist constants 
C, Cl, C 2 > 0 such that 

{UN{t]U)'il:,MUisi{t]U)'ip) < C'exp(ci exp(c2|t|))('i/^, {M PM'^/N + 'Hn) V’) 
\{UN{t-,U)il}, {CN{t) - r]N{t))UN{t-,U)il})\ < (7 exp (ci exp (c2|t|))('(/’, {M + /N PUm) i) 

To prove Theorem 4.2, we use the following lemma, whose proof can be found in [6]. 

Lemma 4.3. Let the fluctuation dynamics U]sr{t-, s) be defined as in (2.15). Then there 
exists a constant C > 0 such that 

W^{t- Q)M‘^ UNit; 0) < C{NU^{t; Q)MUNit- 0) + N[M +1) + {M+ 1)^). 

Proof of Theorem 4-2. From Theorem 3.1, (4.1) and Lemma 4.3, we obtain 

^ {Un{ t; 0 ) 1 ), MUn {t; 0)?/') 

= {lLNit]0)'ijj,[iJ\f, CN{t)]UNit]0)'ijj) 

< Ce^^^\l(N{t;0)f;,AfUNit]0)'il;) + (^^Af(t; 0)V>, {CN{t) - r]N{t))UN{t-,U)ip) 

+ {i), [(AT + 1) + M^/N] i;) 

and 

{CN{t) - r]N{t))L(N{t]0)'ip) 

= {lLN{t;0)ijj,{CN{t) - fjN{t))l(N(t-,0)if) 

< Ce^^^\llN{t;0)'i),MUN{t]0)ip) + {UN{t;0)'i), {CN{t) - r]N{t))UN{t-,U)i) 

+ {if, [(AT + 1) + M^/N] if) 


48 


We conclude that 
d 


— (UNit]0)'ilJ, + 1) + {CN{t) - mit)) UN{t;0)'ilj 


< (l{N{t;0)'ip, Ce^^^^{Af+ 1)+ {CN{t)-VNit)) llN{t;0)'ip^ 

+ (V^, [{M + 1 ) + WViV] V-) 

for appropriate constants C,K (chosen so, that the operator Cexp{K\t\){J\f + 1) + 
{C]y{t) — is non-negative) and D > 0. Gronwall’s Lemma implies that 

^Z^Ar(t;0)V', +1)+ {CN{t)-r]N{t)) 0)V’^ 

< Cexp(ci exp(c 2 |t|)) {ip, [(W -|- 1) +J\f‘^/N + Tiisf] ip) 

which completes the proof of the theorem. □ 


Finally, we control growth of number of particles, energy and their squares w.r.t. 
quadratic fluctuation dynamics. 

Theorem 4.4. Under the same assumptions as in Proposition 2.2, there exist constants 
C, Cl, C 2 >0 such that 

{l{2,Nit;0)ip,MU2,Nit;0)ip) < C'exp(ci exp{c2\t\)){ip, {Af + l)ip) 
{l(2,N{t-,0)ip,J\f‘^U2,N{t]0)ip) < C'exp(ciexp(c 2 |t|))(V', (W-F if ip) (4.4) 

{U 2 ,N{t-, 0 )ip,Cl^Nf)^ 2 ,N{t-, 0 )ip) < Cexp{ciexp{c 2 \t\)){ip, (JC +M + l)‘^ip) 

Remark. The last bound in (4.4) , combined with (4.3), also implies that 
{1^2,N{t', 0)ip,}C'^U2,N{t', 0)ip) < C'exp(ci exp(c 2 |t|))(V', {Ki + Af + l)‘^ip) 

Proof. From Theorem 3.1, we obtain 

-^{h(2,Nit;0)lp, {Af + l)U2,N{t',0)lp) = PU2,N{t',D)lp, [iN,C2,N{t)]U2,N{A, 2l)lp) 

< C'e^l*l(iY2,Ar(t; O)^', {Af + l)Z^2,Ar(t; o)^') 


Gronwall’s lemma implies that 

{l{ 2 ,N{t-, 0 )ip,AfU 2 ,N{t; 0 )ip) < Cexp{ciexp{c 2 \t\)){ip,{Af + l)ip) 

The estimate for Af^ follows analogously. Finally, we compute, using again the 
bounds in Theorem 3.1, 

^ {U2,N{t-, 0)lp, ClN{t)U2,N{t-, 0)lp) 


= {lf2,N{t;0)lp, C2,N{t)C2,N{t) +h.C. U2,N{t;0)lp) 

< {l{2,N{t-,0)lp,Cl^Nit)l^2,N{t;0)lp) -h {ll2,N{t-,0)lp,Cl^N{t)^2,N{t-,0)lp) 

< {h(2,N{t; 0)lp, Cl N{t)U2,N{t-, 0)lp) + Ce^\^\{U2,N{t-, 2P)ip, {Af + l)^Z^2,Af(t; 2i)ip) 

< {U 2 ,N{t\^)ip,if\Nf)'^ 2 ,N{t-,D)ip) -t-C'exp(cexp(c|t|))(V', (Ar +l)^^') 

Gronwall’s lemma and equation (4.3) imply the last bound in (4.4). □ 


49 










With the help of Theorem 4.2 and Theorem 4.4, we can now conclude the proof of 
Proposition 2.2. 

Proof of Proposition 2.2. We have 


0)Cn — e *-^0 0 )^Ar 


dt 


= 2lm(UN{t,0)f,N,{^N{t) — ^2,N{t) —dN{t))e 

< C\(fAN{t,^)iN ,£{t)U2,N{t]^)iN)\ + C'|(Z^Ar(t, 0)^7V, V7vZ^2,Ar(t; 0 )^ 7 v)| 




{hlN{t]t))^N-, {KL +Vn + + '^)UN{t]^)iN) 

+ (^2,Af(i; 0)CAr, {K? + + l)^^2,Ar(t; 0)^Ar) 


where a = min(/3/2, (1 — /3)/2) and where we used the bounds in Theorem 3.1 and 
Lemma 3.10. Applying Proposition 4.1, and then Theorems 4.2 and 4.4, we obtain 

— UN{t] 0),^Ar — e“*-^o 0)^Ar 

< C'exp(ci exp(c 2 |t|))A^“"(^Ar, + /C^ + VAr)^Ar) 

Integrating over time, we conclude the proof of Proposition 2.2. □ 


5 Comparison with the limiting fluctuation dynamics 


Next, we prove Theorem 1.2. Proceeding similarly as in Proposition 2.2, we observe that 




= W{^/Nip^)TN,t K(i;0)^7v -e-'/o^iv(A<i-z^2,oo(i;0)eiv 


From the result of Theorem 1.1, it is enough to compare i^ 2 ,Ar(t; 0)^7v with Z^ 2 ,oo(i; 0)^Ar. 
To this end, we need to compare the two generators C 2 ,N{t) and C 2 ,oo{t) dehned in (1.31) 
and (1.40); we do so in the next four lemmas. 

Lemma 5.1. Recall the definition of and C^^\t) given in (1.33) and after (1.40). 
Under the same assumptions as in Theorem 1.2, with a = min(/3/2, (1 — /3)/2), we have 

(V'l, (45 (0 - '^ 2 ) 

< CN~°‘exp{ci exp(c 2 |t|)) \\{M + l)^'^^'0i|| ll(AA + 1)^'^^V’2|| 

Proof. From (1.33), all contributions to the difference have the form 



Al = J dxa*a^{mx) 

^ 2 = [ a^{jx)a^{mx) 


(5.2) 
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Here m, j S x M^) with ||j ||2 < and 

||m ||2 < exp(ci exp(c 2 |t|)) (5-3) 

Eq. (5.1) follows therefore from the bounds 

\{i)i,Ai'ip2)\ < C'||m||2||(AA+ ||(Ar+ 

and 

KV'l,^2-02)1 < C'||m|| 2 ||j|| 2 ||(AA + ||(Ar+ 1)^/2'02|| 

which can be proven as in Lemma 3.4. 

When we consider the different terms in (1.33), the kernel m is a difference like 
kN,t - h, PN,t - Pt, rM,t - n, VpN,t - '^Pt, ^PN,t - ‘^pt or Ar]sf,t - An. In these cases, 
(5.3) follows from the results of Appendix C. 

Notice that also the contribution 

B = J dx\7,a*{k^)VMk^)- j dxV^a*ik,)VMkx) (5.4) 

can be written in the form Ai in (5.2), with m = — u and 

UN,tix,y) = j dzVxkN,tix, z)VxkN,tix,y), ut{x,y) = j dzVxkt{x,z)V^kt{x,y) 
To prove that (5.3) holds in this case, we estimate 
\\uN,t - utWl = j dxdy\uN,ti^^y) “ utix,y)\‘^ 

— ^ J dxdy / dzVzkN,tiz,y)iVzkN,tix, z) - V 2 /cf(x, z)) (5.5) 

+ C J dxdy J dz{VzkN,t{z,y)-Vzkt{z,y))V Jit{x,z) 

Integrating by parts in the first term, we obtain 

J dxdy J dzVzkN,t{z,y){VzkN,t{x,z) - Vzkt{x,z)) 

< I dxdy / dz AzkN,t{^,y){kN,t{x,z) - fcf(x,2:)) 

< j dxdydzidz2 \Azj^kjs[,t{zi,y)\\Az2kN,t{z2,y)\ 

X \kN,ti^, Zl) - hix, Zi)\ \kN,t{x, Z2) - kt{x, Z 2 )\ 
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By Cauchy-Schwarz inequality, we find 


jdxdy / dzVzkN,t{z,y){VzkN,t{x.,z) - V^A;f(x,z)) 

< / dxdzi\kN,t{x,zi) - kt{x,zi)\^ svip / dy | A^j/cAr^f (zi, 2 /)| sup / dz 2 \i^z 2 kN,t{z 2 ,y)\ 
J z\ J y J 


<Ce^\%kN,t-kt\\l 


In the last inequality, we observed that 
Az^kN,t{^2,y) 

= - - y)iv^{A2 + y)/Af 

N 


(5.6) 


- -^^znA^ 2 - y) AtiA2 + y)/2)A(p^{{z2 + y)/2)) + {Vipff{{z 2 + y)/2)\ 


^^{{z2 + y)/2)A^f{(- ' 

- 2 NVujnA ^2 - y)^f{{z2 + y)/ 2 )V(pf{{z 2 + y)/ 2 ) 
Using the scattering equation for = 1 — we get 
Az:^kNA^2,y) 

]V30 

= A2 - y))fNA^2 - y)i‘ff^{{z2 + y)/Af 

^N, 


- N\N/fNA^2 - y)A\z2 -y\< t){At{A2 + y)/Af 

N 
~2 


- yWAr,£(^2 - y) At{A2 + y)/ 2 )Aif^{{z 2 + y)/ 2 )) + {Vip^f{{z2 + y)/ 2 )] 

- 2NVujnA^ 2 - y)(pfi{z2 + y)/2)Vipf{{z2 + y)/ 2 ) 

With the bounds from Lemma 2.1 and Proposition B.l, we conclude that 

sup [ dy\AkNA^ 2 ,y)\ < 11^2 < 

Z 2 J 

The results of Appendix C, combined with (5.6), imply that 


dxdy 


dzVzkNA^^ y)(yzkN,tix, z) - V^kt{x, z)) 


<CN " exp(ci exp(c 2 |t|)) 


Proceeding analogously to bound the second term on the r.h.s. of (5.5), we obtain 


- uth < Cexp(ci exp(c 2 |t|)) 


_l_ ]y-d/2 


Hence, as claimed, also the term (5.4) can be written as the term Ai in (5.2), with the 
kernel m = ujy — u satisfying (5.3). □ 
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Lemma 5.2. Recall the definition of d'^N given in (1.32) and after (1.40). 

Under the same assumptions as in Theorem 1.2, with a = min(/3/2, (1 — /3)/2), we have 

(V'l, - dp {t)'j fi2) 

< C'7V“"exp(ci exp(c 2 |t|))||(AA + l)^/^'0i|| ||(AA +/C + 1)^/^V’2|| 

(In fact, the kinetic energy operator K. could also be placed on fii.) 

Proof. From (1.32), we observe that all terms in C^pt) have one of the following forms: 


^1,1 = j dxdyN^^V{N^ix - y))(ppx)p{y)a^{jP)a^{jP) 

^1,2 = J dxdyN^^V{N^{x - y))ippx)p{y)a\iP)ay 

= J dxdyN^f^V{N^{x - y))p{x)p{y)alay 

= I dxdyN^^V{Ndx-y))P{x)P{x)J{pydHj2d 
^2,2 = J dxdyN^^V{Nl^{x - y))p{x)p{x)a^{py)ay 

= j dxdyN^I^V{Ndx - y))<Pt {x)p{x)a*yay 


(5.7) 


(or possibly, the form of the adjoint of A 12 or A^o)- Here ip^ denotes the solution 
of the A^-dependent nonlinear Schrodinger equation (2.10), and 3^)2 ^ L^(M^ x M^) 
being either the operator sinh^^^ or pw.t = cosh^j^ ^ —1. (In fact, some of the factors 
should be replaced by pf, but this does not affect our analysis). 

To estimate the difference dp{t) — dp{f), we need to compare the terms in (5.7) 
with their formal limits: 


7li,i = 
^ 1,2 = 
^ 1,3 = 
^ 2,1 = 
^ 2,2 = 
^ 2,3 = 


5 o J dxdy 6 {x-y)ptix)ptiy)a^{ji,x)d{j 2 ,y) 

bo J dxdydix - y)pt{x)pt{y)a^{ji,x)ay 

60 J dxdy 5{x - y)pt{x)pt{y)alay 

60 J dxdy6{x - y)pt{x)pt{x)a^{ji^y)a^{j 2 ,y) 
bo J dxdyd{x - y)pt{x)pt{x)a^{ji^y)ay 
bo J dxdy6{x - y)pt{x)pt 


ot{x)a*a. 


where pt is the solution of the limiting nonlinear Schrodinger equation (1.26) and ji, j 2 £ 
L2(R3 X m 3) are either sinh(A:t) or pt = cosh^^ —1, with kt given by (1.37). Note that, 
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from the results of Appendix C, we always have 

\\ji - 3ih < C'iV"“exp(ciexp(c 2 |t|)) (5.8) 

for i = 1,2. The lemma will follow if we can prove that 

< C'lV“"exp(ci exp(c 2 |t|))||(Ar + ||(AA + /C + 


for all i = 1,2 and j = 1,2, 3. 

We start comparing with Aip. To this end, we observe that 

= j dxdyN‘^<^V{N^{x - y))<ft {x)^t {y)aK3i,x)aK32,y “ hy) 

+ I dxdyN^^V{Nl^{x-y))ipf{x)ip^{y)J{j^^^-ji,x)a^{j2,y) 

+ j dxdy N^^V{N^{x - y))-ho6{x - y) {x)ipf {y)J{ji^x)a^ij 2 

+ j dxdybo6{x-y)(^y?^{x)y?^{y)-iptix)(ptiy)^J{ji^x)a\j2,y) 

The contribution from the first term can be bounded by 

f dxdyN^l‘V(N^(x - faXV'i, - i2,j)V'2) 

< C\\y,”f„, ||jf lb 11*" - hh ll(V+ 

< CA^““exp(ci exp(c 2 |t|))||(AA + ll(A^ + 1)^^^^2|| 


(5.9) 


(5.10) 


Also the contribution from the second term can be bounded similarly. The fourth term 
on the r.h.s. of (5.9), on the other hand, is bounded by 


/ 


dx {'ipi,a^{ji,x)aHj 2 ,x)'ip 2 ) 

< C'dlv^f Iloo + ||<y2t||oo)||‘/?f - V^tlloo j dx\\J{ji^x)'ipi\\\\aKj2,x)'ilJ2\\ 

< CA“"exp(ciexp(c 2 |t|)) \\{Af + lf^^ipi\\\\{Af + l)^^‘^'ip 2 \\ 

Finally, we have to bound the contribution arising from the third term on the r.h.s. of 
(5.9). We have 


dxdy {N^{x - y)) - boS{x - y))(pf ix)(pf {y){'ipi,J{ji^^)J{j 2 ,y)'ip 2 ) 

= j dxdyV{y) (x)(/?f (x + y/A^^)('0i, at*(ji,a;)a>*(j2,,,+j^/7V/3)V'2) 

- (a;)^('0i,«‘*(ji,x)a**(i2,x)'02) 
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and thus 

j dxdy {N^^V{N^{x - y)) - bo6{x - y))^p^{x)ipf {y){'tpi,J{ji^x)a^{j 2 ,y)^ 2 ) 

< J dxdyV{y)\ip^{x)\\ip^{x + y/N^) - ip^{x)\\\J{ji^^)i;i\\\\J{j 2 ^^+y/Nf^)'ip 2 \\ (5.12) 
+ j dxdyV{y)\ip^ix)\‘^\\J{ji^^)'ilJi\\\\J{j2^^^y/Nfi - j2,x)ilJ2\\ 

The first term can be easily controlled with Proposition B.l. We find 

J dxdyV {y)\ip^ {x)\\ip^ (x + y/N^^) - </?f (x)|||a“(ji,x)V’i||||a“(j2,x+y/7V/3)V’2|| 

< CN~^ exp{ci exp(c 2 |t|))||(AA + l)^/^i/>i||||(Ar + 1)^/^V'2|| 

As for the second term, we estimate it by 

J dxdyl/(y)|(/?f (3:)p||a**(ji,a;)V'i||||a‘‘(j2,x+y/7V/3 -j2,x)V’2|| 

< II (AA + II II (AA + II / dxdy 

X f dXV{y)\y\\ipf{x)\'^\\ji^^\\2\\V:,j2^^^XyN-42 

Jo 

<CA^-^e^l*l||il||2||V2j2||2||(AA+l)l/2^l||||(AA+ 1)^2^211 

We are interested in j 2 = sinh^j and j 2 = Pt = cosh^j —1. In both cases we have 
IIV2J2II2 < (in the case j2 = Pt, we actually have the better bound ||V2j2||2 < C*; 

see Appendix C). Hence, we conclude that 

J dxdyV{y)\p^ -i 2 ,x)V' 2 || 

< CiV-^/2gi^|t| II ^ 1)1/2^^ II II ^ 1)1/2^2 II 

and therefore, inserting in (5.12), that 

j dxdy{N‘^^V{N^{x - y)) - bo5{x - y))Pt {x)ipt {y){i}i,a^{ji,x)a^{j 2 ,y)'^ 2 ) 

< C'A^“^/2gxp(ciexp(c 2 |t|))||(AA + 1)^/2^^|| ||(Ar+ lY/‘^'ip 2 \\ . 

Together with (5.10) and (5.11), we obtain from (5.9) that 

|(V’i,(A(^i -^i,i)V’ 2 )| < C'A^"“exp(ciexp(c 2 |t|))||(AA + l)^/^V’i|| ll(AA+ 1)^/^V'2|| 
Similarly, we can show that 

|(V’i,(A (^2 -^i, 2 )V’ 2 )| < C'Af"“exp(ciexp(c 2 |t|))||(AA + l)^/^V’i|| ll(AA+ 1)^'^^V'2|| 
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In fact, the contribution 

j dxdy (^N^<^V{Nf^{x - y)) - boS{x - y)^ {x)(p^{y){ipi,J{ji^x)ayip 2 ) 

< J dxdyV{x) pf{y)pt{y + xN~f^){'4}i,a\ji 

^y^xN~^ )ay'ljj2) 

< [ dxdyV {x)\p^ {y)\ pt{y + xN-^) - (y) \\a^iji^y+^N-f>)'4’i\\hyip2\\ 


+ 


J dxdyV{x)\p^{y)\‘^\\J{ji^y^^]^-i 3 - ji,y)ipi\\\\ay'ilj 2 \\ 


can be dealt with as the third term on the r.h.s. of (5.9). 

Next, we consider the difference — ^ 1 ^ 3 . We write 

- ^i,3)V’2) 

< j dxdyN^^V{N^{x-y))\pf{x)-pt{x)\\pf^{y)\\\a^'ipi\\\\ay'ip2\\ 

+ J dxdyN^^V{N^{x - y))\ptix)\\p^ (y) - pt{y)\\\axM\hy'^2\\ 

+ j dxdyN^^V{N^{x - y)) [pt{x)pt{y){axi}i-,ay'il} 2 ) - Pt{xf {axil}i,a^'4;2)] 

= I + II + III 

(5.13) 

We can control the first term on the r.h.s. by 

I < \\Pt - ¥^i||oo||¥^t||oo||(Ar+ l)^^^^i||||(Ar+ 1)^^^V'2|| 

< C'IV“"exp(ciexp(c 2 |t|))||(AA + l)^/^V’i||||(-^ + 1)^'^^V'2|| 

The second term on the r.h.s. of (5.13) can be controlled analogously. As for the third 
term, with a change of variables we can write 


III < 


J dxdyV{y)\ptix)\\pt{x + yN ^) - ptix)\\\ax'4’i\\\\a^+yj^-p'ilJ2\\ 

J dxdyViy)\pt{x)f\{ax^l^i 

1 i^x+yN~P )V’2)| 


< C'e^l'llV-^||(AA + 1)1/2 ^i||||(AA + 1)^/2V;2 
+ J dxdy J dXV{y)\y\N~f^\pt{x)\'^\\a, 


(5.14) 


'^xax+\yN-f^'^ 2 \\ 


< Ce^l*llV-^||(Ar+ l)i/2^i||||(/C +AA + 1)^/V2| 

The bounds for the terms 

|(V’l,(^^i -^ 2 ,i)V’ 2 )| 


(5.15) 
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for i = 1,2,3 can be obtained similarly. In this case, since both field operators depend on 
the same integration variable y, the difference between {N^{x — y)) and bQ6{x — y) 

can be controlled through the difference between [x) and [x + yN~/^). Smallness 
here follows from the regularity of , there is no need to use the kinetic energy of '01 
and 02 (in contrast with the term — ^ 1 , 3 ). To illustrate this point, let us briefly 
consider the term 


( 01 ; (^2,3“^2,3)02 

< 


J dxdy{N^^V{N^{x-y)) - 6 o 5 (x - y)) yjf (x)^(ay0i, aj^02) 

+ bo f dx (x)p - |v9i(x)l^| ||a3;0i||||a3;02|| 


(5.16) 


The second term can be bounded by 

^0 ||co + llv^tlloo] llv^t ~ y^illo 

< C'A^““exp(ci exp(c2|t|))||Ar^/^0i||||Ar^/^02 


1/2 / . ^ 1/2 

dx ||aa;0i|p ) { dx ||aa; 02 |p 


As for the first term on the r.h.s. of (5.16), it can be estimated by 
j dxdyV{x) ip^ {y + xN-f^f - if^ {yf ||aj/ 0 i|| ||ay02|| 

j dxdy J (iAl/(x)|x||V</?f (y + AxA^“''^)|||aj/0i||||ay02|| 

From the last two estimates, we conclude that 

1(01, (A^3 - A 2 , 3 ) 02)1 < C'A^""exp(ciexp(c2|t|))||AA^/^0i||||AA^/202|| 

Analogously, we can also control (5.15), with i = 1,2. 


□ 


Lemma 5.3. Recall the definitions (1.29) and (1.39) of T]\f^t o,nd Tt. Under the same 
assumptions as in Theorem 1.2, with a = min(/3/2, (1 — (3)/2), we have 

1 (01, {{idtTf,^fiTM,t - {idtTm) 02)1 

< Cexp(ci exp(c 2 |t|))A^"“||(Ar+ l)^/^0i||||(AA + 1 )^/^ 02 || 

Proof. We use the notation 

= ^ j dxdy{kN,tix,y)ala*y-kN,t{.x,y)axay 

and ^ 

B{t) = T, j dxdy[kt{x,y)ala* -kt{x,y)axa^ 
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Then we have Tw,* = exp{B]\f{t)) and Tt = exp(i?(t)). We will make use of the repre¬ 
sentations 


with 


adsj^(t)(-B7v(i)) = 2 j dxdy {f^^i{x,y)a*^a*y + fj^^2{x,y)a^ay) 

adi3(t)(^(i)) = ^ j dxdy {fn,iix,y)alal +fn,2ix,y)a^ay) 

where and fn^i, for i = 1,2 and n € N satisfy the bounds from Lemma 3.12 (once 
with kjsf^t and once with the limiting kernel kt). We find 

KV^i, - [idtT:]Tt)^P2)\ 

1 


< 


^ I lu K^i’ (adB^(t)(57v(t)) - ad^(,)(5(t)))V^2>| 


< 


n>0 

E 

n>0 

+ IIV'i 


1 


(n + 1)! 


(Il/h - fn.lh + \\&- /„,2||2)|I W + l)'/LlllllW + 1)‘''L2| 


EtoW I |/ 2 Li, 2 (*.*) -/2n-l,2(a:,a:)| 

n>l '> 


(5.17) 


Next, we claim that 

fnA\2<[2\\kt\\2T\\kN,t-h\\2 


n—1 


+ 2\\kN,t - kt\\2\\kN,t\\ ^ {A\kN,t\\y {A\kt\\) 

j=0 


n-l-j 


(5.18) 


and that 


J dx\fyi{x,x) - fn,i{x,x)\ < {2\\kt\\2)'^\\kN,t - hh 


n—1 


+ 2\\kN^t - ^( 2 ||t„||)'( 2 ||t,||) 

2=0 


n-l-j 


(5.19) 


for all n G N (with the convention that the sum disappears if n = 0). These bounds, 
together with the estimates in the appendix C, complete the proof of the lemma. 

It remains to show (5.18) and (5.19). We proceed by induction over n G N. For 
re = 0, we have 

y) = kN,t{x, y), fi^ 2 {x, v) = kN,t{x, y) 
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and analogously for /o,i and /o, 2 ; the bounds (5.18) and (5.19) are clearly satisfied. Let 
us now assume that (5.18) and (5.19) hold for some n G N. We prove they also hold for 
(n + 1). We consider first the case n even. By Lemma 3.12 we have 

/„+i,i(x,z) = -^ J dy(ktix,y){fn, 2 iz,y) + fn, 2 iy,z)) 

+ h{y, z) {fn, 2 ix, y) + /n, 2 (y, x))j 
fn+i, 2 (x,z) = -^ j dy(ktiy,z){fn,i{x,y) + fn,iiy,x)) 

+ h{x,y){fn,iiz,y) + fn,i{y,z))'^ , 

and similarly for i{x, z) and f^+i 2(^1 ^)- Therefore 

211/2 

\\fn+l,l - /n+l,l||2 < 



1 

+ 2 


+ 


J dxdz j dy{kN,t{y^z)fn,2ix,y) - kt{y,z)fn,2{x,y)) 
j dy{kN,t{y, z)fn^2iy, x) - kt{y, z)fn,2{y, x)) 


dxdz 


11/2 


1/2 


11/2 


By the Cauchy-Schwarz inequality, we find 

ll/n^+ 1,1 - /n+l,l||2 < 2 ( 11 %,t - kth\\f^,2h + ll^ilbll/^ - /n,2||2) 

Using the induction assumption (5.18) and the bound ||/^j ||2 < (2||A:Ar,f || 2 )"'||^Af,t ||2 from 
Lemma 3.12, we obtain (5.18) with n replaced by (n -|- 1). Moreover, we notice that 

J dx \f^^i^^i{x,x) - fn+l,l{x,x)\ 

= ^jdx J dy{kN,tix,y) + kNAy^^))ifn,2ix^y)+fn,2iy>^)) 

- J dy{kt{x,y) + kt{y,x)){fn,2A,y) + /n,2(2/,a;)) 

< jdxdy\kN,tA^y)-ktix,y)\{\fA,2A,y)\ + \fn,2iy^^)\) 

+ 2 j dxdy\kt{x,y)\\fAAx,y) - fn,2A,y)\ 

< 2 \\kN,t - hUfA,2h + nfn,2 - fn,2h\M2 
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Again, the induction assumption, combined with the bounds of Lemma 3.12 imply (5.19), 
with n replaced by (n + 1). The bounds for i = 2 and for n odd can be proven similarly. 

□ 

Finally, we compare the last two terms on the r.h.s. of (1-31) with the corresponding 
terms on the r.h.s. of (1.40). 

Lemma 5.4. Under the same assumptions as in Theorem 1.2, with a = min(/3/2, (1 — 
/3)/2), we let 

= J dxdy Nujn,^{x - y)(j)N{ix + y)/2)a*a* 

D 2 = N>^n,£ j dxdyl{\x - y\ < l)Tt + y)/2)^<Oy 

and 

Di = y dxdy {x - y)(t>{{x + y)/2)a*al 

-^2 = j dxdyl{\x -y\< i)Mix + 2 /)/ 2 )^a*a* 

with 


4>Nix) = (p^{x)A(pf^{x)+ \V(pf^{x)\‘^ and 4){x) = <pt{x)A(pt{x) + \Viptix)\‘^ 
Then we have 

KV'i, (A,AT - A)^ 2 )| < C'iV““exp(ciexp(c 2 |f|))||(AA + 1)^/^AIIII(A^ + 1)^'^^V'2|| 


for i = 1,2. 
Proof. We have 


Df - Di = J dxdy {Nu}e,N{x - y) - - y))(l)N{x + y) 

+ J dx dy {x - y) {4 >n{x + y) - (j){x + y)) alal 


* * 


(5.20) 


To estimate \ {D^ — Zli)^)!) we observe that, for any Hilbert-Schmidt operator A 

on L^(M^), with the integral kernel A{x,y), we have (with the usual notation Ax{y) = 

A{x,y)) 

J dxdyA{x,y){ifi,ala*yil)2) = J dx(aa;A, «*(^a;)^2) 

J (ix||aa;AIIPx||2||(AA +l)^/^i/>2|| 


< 

< 


hs||(AA + 1)V2^i||||(A+1)V2^2 
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Hence, we find 


<||(AA + 1)V2^i||||(AA + 1)V2^2|| 

X J dxdy\NujN,e{x - y) - - y)f\(l)N{{x+ y)/ 2 )\'^ 

+ j (ixdya;f^™P(x-y)^|(/)Ar((x + y)/2)-(/)((x + y)/2)| 

< C'iV""exp(ci exp(c 2 |t|))||(AA + IIll(AA + 1)^/^V'2|| 

To control the difference — D 2 , we write 

{D^ - D 2 ) = {N\,,n - Xi) J dxdyli\x-y\<i) ((x + y)/2))" 

+ Xi J dxdyl{\x-y\ < £) ((x + y)/ 2 ) )^ - + 

where we defined = Sbo/Sni^. Hence, with (5.21), we find 

- 02 )^ 2 ) 

<C||(AA + 1)VVi||||(AA + 1)V2^2|| 


11/2 


NXni — 


36o 


SttP 


+ C'e^l*l||(/?f-¥5t||2 


< CN "exp(ci exp(c 2 |t|)) ||(A/ + l)^/^V’i|| ll(AA + 1)^/^V’2| 


* * 


□ 


We are now ready to show Theorem 1.2. 

Proof of Theorem 1.2. We write 

\e-^'^^^W{VN^)TNfliN - e-^foVN{s)ds W(ViV(^f )r^,tZ^2,oo(/)C7v| 

< ||Z^7v(/; 0 )^Ar — e“*-^o 0 )^Ar|| 

< ||Z^iv(t; 0)^Ar — e“*-^o 0)^Ar|| + \\U2,N{t; 0 )Cn — ^2,oo(/; 0)^Ar|| 

The first term can be estimated with Proposition 2.2. To bound the second term, we 
apply Lemmas 5.1, 5.2, 5.3, 5.4. We find 

^||^2,Af(L 0)'^Af -Z^2 ,oo(L0)^1v|P 

< C\(U2,Nit]0)f,N, {C.2,N{t) — T2,cx)(i))^2,Af(/;0)^Af)| 

< C'Ai““exp(ci exp(c 2 |t|)) 

X [(iY2,Ar(t; 0)^Ar, (A/ + /C + l)/Y2,Ar(t; 0).^Ar) + (Z^2,oo(/; 0)^7V, (A/ + 1)Z^2 ,cx)(L 0)?Af)] 

< C'Af““exp(ci exp(c 2 |t|))(CAr, (A/ + /C + l)CAr) 
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Theorem 1.2 follows integrating this bound over time. In the last step, we used Theorem 
4.4 and the fact that, exactly as in Theorem 4.4, we can also control the growth of the 
expectation of J\f with respect to the limiting fluctuation dynamics Z^ 2 ,oo, he. there exist 
constants C,ci,C 2 >0 such that 

< C'exp(ciexp(c2|t|))(CAf, (AA+ l)^Ar) 


□ 


A Properties of the Scattering Function 


In this appendix, we give a proof of Lemma 2.1, containing basic properties of the ground 
state /tv,^ of the Neumann problem 




In/ = ^N,efN,e 


on the sphere |x| < i. 


Proof of Lemma 2.1. To prove an upper bound on Xn/-, we use a constant trial function. 
With the notation 

\) = -A+^N^^V{Nl^x), 


2N 


we find 


Xnj < 


( 1 , 1 , 1 ) _ 1 _ 3bo 


(1,1) 2N 




\x\<i 


dx 


SttNP 


for i ^ N~^. Before proving the lower bound in (2.6), we prove parts ii) and iii) of the 
Lemma. 

Part ii) can be proven as in [16, Lemma A.l]; we skip the details. As for part iii), 
we observe that, from (2.5), we can write 


0 JN,i{x) = C J dy ^y{Ady)fN^i{y) - XN,ifN,i) 


< 


c 

iv 


/ dy^^iv3V(iv/'y)/w.Ky) 

J \x-y\ 


(A.l) 


The Hardy-Littlewood-Sobolev inequality implies that 


C 


(jxj +N d)u)]si^i{x) < — dy 


+ N-d 


< 


1 + 


\x - y\ 
C 


N^dv{N^y) 


N 
C 
N 

<h[||r||i + ||(|.| + i)F(.)M 


f dy-2—N^Sv(Nl>y){N^\y\ + 1) 

J \x-y\ 
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Furthermore, taking the gradient of (A.l), we can also estimate 


\{N + \x\‘^)Vu}N,e{x)\ 

^Cjdy (+ 1 + (jVW-> V(Nl>y} + A„,,l(| 9 | < (j) fy/y) 

< C(||iV^^-V(iV^-)lli + + 1)||3) 

+ CA„,,(iV-"'>||l(|!/| < «)||3 + ||1(|!/| < <)||, + |||</|U(|!,| < «)||3) 

<^{l|r||3 + l|r||, + |||.|M.)ll3 + i) 

for any i ^ N~f^. 

Finally, we have to prove a lower bound for the eigenvalue Xn/- To this end, we use 
the estimates 0 < iON^£{x) < C/{N\x\), established in part ii) and iii). We have 


XN,e = 


hfN,e) 

{Jn/, fN,e) 


> 


fi 


\x\<i 


1 + ujj^^ix )) dx 


2N 


|x|<€ 


N^f^V{Nf^x) (1 - 2ujn/{x)) 


- 4^/3 I « 2iV 
3''^ ^ IW 


c 


bo-j^\\V{.)/\.\h 


> 


36o 

8 ttNP 


1 - 


C 


iVl-/3 


for all £ ^ and 0 < /3 < 1. 


□ 


B Properties of Nonlinear Schrodinger Equations 

For a given initial data (p £ we define as the solution of the modified time 

dependent nonlinear Schrodinger equation 

idtPtix) = -Apl^(x) + (JV^^V(N^-)fN,£(-) * l‘Pt^l^)(x)pf^(x) (B.l) 

where fN,e is the solution of the scattering equation (2.5). Moreover, we denote by pt 
the solution of the limiting equation 


idtpt = -Apt + bo\pt\‘^pt 


(B.2) 


where bg = f dxV(x). 

We need some standard facts concerning the well-posedness and the propagation 
of higher Sobolev regularity for these equations. Moreover, we need to estimate their 
difference in the limit of large N. This is the content of the next proposition. 
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Proposition B.l. Let V £ n (1 + \x\^)dx) be non-negative and spherically 

symetric. Let (/? G with ||</?||2 = 1- 

i) There exist unique global solutions and (p, in (^(M; of (B.l) and, respec¬ 

tively, of (B.2) with initial data ip. The solutions are such that \\p>t\\2 = IIII2 = 
||(/?||2 = 1 and 

llv^ilUi) \\Tt Wm < C 

for a constant C > 0 and all t gW. 

a) Under the additional assumption that ip G for an integer n G N, then 

ipt,T^ G for all t £ M. and there exist constants C > 0 (depending on 

||(/ 9 ||/fn and on n) and K > 0 (depending only on Ilf/Jlli^i and on n) such that 

llv^tll-ff"! Hi/" < 


for all t G M. 

Hi) Let ip £ Then there exists C > 0 (depending on ||(/ 5 ||i^ 4 j and K > 0 

(depending only on ||(/ 9 ||j:^i j such that 

llv^tlli^a, \\ipt\\2 < 

iv) Let ip £ Then there exist constants C,ci >0 (depending on ||(/?||H 2 j and 

C 2 > 0 (depending only on ||</^||hi j such that 

\\ipt - ipt \\2 < CiV”'’'exp(ci exp(c 2 |t|)) with 7 = min(^, 1 - (3). 

v) Let ip £ H^(R^). Then there exist constants C, ci >0 (depending on ||<^||H4j and 
C 2 > 0 (depending only on ||<^||hi j such that 

\\Tt - \\Tt - T^th < CiV”'^ exp(ci exp(c 2 |t|)) with 7 = min(/3,1-/3) 


Proof. The proof of i)-iii) is quite standard and can be found, for example, in [ 6 , Prop. 
3.1]. Also the proof of the bound 

\\ipt - ipf \\2 < C'/V“'’'exp(ci exp(c 2 |t|)), with 7 = min(/3,1-/3) 
can be found in [ 6 , Prop. 3.1], up to the observation that 


I dxN^^V{Nf^x)fNM - 6c 




as follows from Lemma 2.1, writing = 1 — To prove v), we observe that 

ipt-ipf = i f ds [{U * \ip^\^)ipf - bo\ips\‘^Ts] 

Jo 
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with Un{x) = N^I^V{N^x)f]\f^i{x). We estimate 


II //2 


^ V"vf V"2^f)V“3y,f-6oV"VsV"2^,V"Vs 

I 1,1 Ml iMn-'O 


|ch|+|«2|+|«3|<2 


< 


^ / ds ||(C/7v*(V“Vf-V“Vs)V“2^f)V“Vf|| 


|ai| + |o2|+|a3|<2 ' 


|«l| + |«2|+|a3|<2 ^ 

|ch|+|« 2 |+|o 3|<2 

M 1 , M 1 , M ,^„JO 


|ch|+|« 2 |+|o 3|<2 

Hence, we obtain 


- ‘Pt Il-ff2 


< 


IlH^lli rhs[||V“ 3 <^f||^||V“Vflloo||V“Vf-V“V .||2 

1,1 Jo L 


|ai| + |a2| + |a3|<2 


+ 


E 


+ l|V“V.||oo||V“Vf l|oo||V“Vf - V“Vs||2 
+ ||V“Vs||oo||V“Vs||oo||V“Vf - V“V.||2 

2ds \\{U * V“VsV“2^jV“V, - 


|«i|+|«2|+|«3|<2 ' 


and using the propagation of regularity from part ii) of the proposition, we conclude 
that 

<C re^N||(^,-¥pf||H 2 +CiV- 7 e^N 
Jo 

By Gronwall, we find 

\Wt - 11^2 < C'iV“'^exp(ci exp(c 2 |t|)) 

The estimate for — ^^\\2 follows with the help of the equations for ipt and (pf. □ 

C Properties of the operator kN,t 

We define the integral kernel 

kN,t{x, y) = -Nujn,i{x - y){ipf{{x + y)/2)f (C.l) 
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where = 1 — ^^n,i is the solution of the scattering equation (2.5). Note that is 
the integral kernel of a Hilbert-Schmidt operator on In fact, using part iii) of 

Lemma 2.1, we hnd 


\\kN,t\\ns = j dxdy\kN,tix,y)\‘^ = J dxdyN‘^u!%^^{x - y)\(p^{{x + y)/ 2 )f 
<C I dxdy + y)/2)|4 < ci 


(C.2) 


In the next hve lemmas, we collect some properties of the kernel kjsf^t and of the 
operators SN,t = sinhfcj^_^, pN,t = coshkj^ _^ -1, rN,t = SN,t - kN,t- 

Lemma C.l. Let V be as in Lemma 2.1 and k]\f,t be defined as in (C.l), with ip = 
G (so that for all t G M, by Proposition B.l). Then 


{i) , ||PAf,t ||2 , ||■SA ^,^||2 , lkAf,t ||2 < C 

{a) sup \\kN,t{x, oib , sup \\pN,t{x, oib < , 

a;GR3 xGR^ 

sup ||sAf,t(a:, Olb, sup \\rN,tix, ■)\\2 < 

xGR® xGR^ 

(m) sup|rAr,t(x; 2 /)|, sup |pAf,t(a:, y)l < 

x^y x,y 


Proof. We start with the property (i). The bound for kjsf^t is proven in (C.2) above. The 
bounds for ||p 7 v,t|| 2 ) and ||sAr,t ||2 follow from 

\\K 1 K 2 h < l|iiil||2||i^2||2 

for any two kernels Ki, K 2 G L^(M^ x M^) (equivalently, ||iLiiL 2 ||HS < ||-fiii||Hs||-fii 2 ||HS 
for any Hilbert-Schmidt operators Ki,K 2 ). 

Next, we show (ii). We observe that 


\\kNAx^-)\\2 = J dy\kN,t{x,y)\'^ 

<ldy \v\p^{y)\^f < CWp^WUlpfWlo < 


using Hardy’s inequality. To prove the bound for sup^, \\pN,t{x, .)|| 2 , we use 


\\{kN,tkN,tTA^-)\\2= j dy j dzkN,tA^A[{kN,tkN,tT ^^N,t]A^y) 


< j dydzdwlkNA^^Al"^ \[ikN,tkN,tT ^kN,t] i'>^^y)\^ 

< ||^Af,t(a;, Olli \\ikN,tkN,tT~^^N,t\\2 < \\kN,tA^-)\\2 \\kN,i 


||4n-2 

II2 
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This gives 


sup \\pN,t{^^)h < sup \\kN,t{x, Olb < 

X X 

as claimed. The bounds for sup^. \\rN,t{x, ■)\\2 and sup^. ||sAr^t(a;, .)||2 can be proven simi¬ 
larly. 

Finally, we show hi). From the definition of ttv,*, we find 


rN,t{x,y) = ^ 


^ (2n + l)! 


{kN,tkN,tTkN,tix,y) 


If _ _ 

'^j2riTTV. / [{kN,tkN,tT~^kN,t] {z,w)kN,t{w,y) 

n=l ^ ' 


and therefore 


Iff 

\rN,t{x,y)\ < ——pT / dwdz\kN,t{x,z)\^\kN,t{w,y){^ 

^ {2n + l}l IJ 


1/2 


-|l/2 


< 


yi_ 


dwdz\{kN,tkN,tT kN,tiz,w)\ 


\\kN4f'-^\\kNA^,mkNA-,y)h<Ce^^^^ 


for every x,y ^ The bound for pjyt can be proven analogously. 


□ 


Lemma C.2. Let V be as in Lemma 2.1 and k^^t be defined as in (C.l), with p = 
(/?^o S (so that for all t € M, by Proposition B.l). Then 

(i) ||ViA:,v,t||2, \\V 2 kNA 2 < 

(a) ||ViPAr,i||2, ||V2PAr,i||2, ||ViTAr^tII 2 , ||V2rAr,i||2 < C 

(Hi) ||AirAr_t||2, ||A2rjv,t||2, \\A 1 p 1 sfA 2 , ||A2PAr,t||2 < Ce^^^^ 

Proof. From the definition of k^^t, we find 

Vi/cAr,i(a;, y) = - NVuinA^ “ y){Tt iA + y)/Af 

- NojnA^ - yA^iA + y)A)^Tt{A + y)A) 

Hence 

J dxdy\VikN,tA^y)\‘^ 

+ C J ^|^-^^3^|V(/?f((x + y)/2)|2|(^f((x + y/2)|2 
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where we used the bounds from Lemma C.l. 

To bound the ||VipTv.tlb, ||VirAr^i|| 2 , we need an estimate for ||Vi(A:7v,tfcAr,f)||2- To 
this end, we notice that 


sup / dy\S/xkN,t{x,y)\ 

X J 

I' 


< Csup dy 


\x — yp 

\v>^({x + y)/2)p 

\x — yb 

by Hardy’s inequality. Analogously, we hnd 


b«((l + !/)/2)|2 |y,"((i + y)/2)||Vy^((i + g)/2)| 


\^-y\ 


< C sup / dy 
X J 


+ |V¥.f((x + y)/2)b 


< C 


sup / dx\VxkN,t{x,y)\ < C 
y J 


Hence, we obtain 

j dxdy\VikN,tkN,ti^^y)\‘^ 

< j dxdydzidz 2 \S/xkN,t (a:, 2 : 1 )11 kN,t ( 2 : 1 , y) 11V^, kN,t (x, Z 2 )11( 2 ^ 2 , y)| 

< j dxdydzidz2\VxkN,ti^^ ^i)\\'^xkN,t{x, Z2)\\kN,tiz2,y)\‘^ 


< 11 %, 


t\\2 


sup / dy|V 3 ,%,i(x,y)| 


sup / dx\VxkN,tix, 


< C 


and therefore 


(C.3) 


Now, we are ready to bound ||VipAr_tlb- From the definition, we have 

cxo ^ 
n=l '' 

\\ViPN,th < ||V(%,t%,t)|b < C 

Similarly, one can show the bounds for V2PN,t-, ^i'rN,t and V2'<’N,t- 
As for the estimates involving second derivatives, we claim that 

\\^l{kN,tkN,t)\\2 < 

In fact, 

= - N j dz{/:\u:N,i){x - z){y)f {{x + z)/2)fkN^t{z,y) 

-2N J dz{VujN,i){x - z) [Va;(¥?f ((x + z)/2))2] kN,t{z,y) 
-N dzujN/ix-z) [Aa;(<y9f ((x + 2;)/2))^] kN,t{z,y) 


(C.4) 


(C.5) 
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and therefore 


/^xkN,tkN,t{x,y) = - ^ j dz{Vu}N/){x - {{x + z)/2)fVzkis[,t{z,y) 

- 3A^y dz(yujN,i)ix - z) [V3;(v7f ((x + z)/2)f] kN,t{z,y) 

- j dzujN/ix - z) [Aa;(<y9f ((x + z)/2)f] kN,t{z,y) 

= I + II + III 


where we used integration by parts. The L^-norm of first term can be bounded by 

l(|x — Zl\ < ■^)l(|x — 221 < i) 


IIII 2 ^ J dxdzidz2 


\X — Zl\^\X — Z2\ 

X {{x + Zi)/2)\‘^\(ff^{{x + Z2)/2)\‘^ f dy\Vzj^kN,tizi,y)\\V^ 2 kN,t{^ 2 ,y)\ 

(C.7) 


Notice that 

j dy\y^,kN,t{zi,y)\\^^2kN,ti^,y)\ 


/■ 


< / dy 


+ \Vip^{{zi + y)/2)\-^ 




- yp 


+ 1 


|^f((z 2 + y)/ 2 )h 

lz2 - 2/P 


+ |V<^f((z 2 + 2 /)/ 2 )|^ 


Zl - Z2 


Inserting this bound on the r.h.s. of (C.7) and shifting the integration variables, we 
conclude that 


III 2 < 


< Ce^l ‘1 


dzidz2 


l(|zi| < ■^) 1 (| 22 | < 

kiPk 2 P 


1 

\zi - Z2 


+ 1 


X sup f + zi/2 )Pi+ 2 : 2 / 2 )p 

21,-22 J 


with the Hardy-Littelwood-Sobolev inequality. The second term on the r.h.s. of (C.6) 
can be controlled by 

||II||i < C'||A:7v,t||i sup [ dx - - [\V{{x + zi)/2)\^ + |vjf ((x + 2 ;i)/ 2 )p] 

Zl J \X ~ Zl\ 

X sup [ dZ2-. - 12 ((x + Z2)/2)\^ + ((x + Z2)/2)\^] 

X J \x — Z 2 \ 
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Finally, the third term on the r.h.s. of (C. 6 ) is bounded by 

1 


7 


nilla < c dxdydzidz2-. 


-\kN,t{zi,y)\\kN,t{z2,y)\ 


X — Zl\\X — Z2\ 

X + Zi/2))||,,f{(x + i,i)/2)| + IVvf((x + Xl)/2)P] 

X [|A»>f ((x + X 2 / 2 ))|kf'{(x + 22 )/ 2 )| + IV,,''((x + 22 )/ 2 )| 2 ] 

|»)f((x + z,)/2)|2 


< C'll^wxllz sup / dx |A(pf((x + 2i/2})|S + 


\X — Zl\ 


+ 


|V(^f((x + zi)/2)| 


2 -, 


\X — Zl\ 




X sup / dZ2 

X 


((® + 2 : 2 / 2 ))!^ + 


|(/?f ((x + 22)/2)P 


F - 2 : 2 ! 


+ 


|V(^f((x + 22)/2)| 


2 -, 


F - Z2\ 


< Ce^l‘1 


This concludes the proof of (C.4). From (C.4), we immediately obtain 

||Aip^,t||2 < Ce\M^\\^{kN,tkN,t)h < 

and 

||Airjv,t||2 < C\\kN,the\\^^'^HHkN,tkN,t)h < 
as claimed. □ 

Lemma C.3. Let V be as in Lemma 2.1 and kM,t be defined as in (C.l), with (p = 
G (so that for all t G M, by Proposition B.l). Then 

(i) \\kN,t\\2 < 

(a) \\kN,th < 

(Hi) \\PN,t\\2, pAr,t||2 < 

Proof. The bounds (i) and (ii) follow from 

\\kN,t\\2 = ^ j dxdy\NujpM{x-y)\^\p^{{x + y)/2)\^\fi(^{{x + y)/2)\^ 

< C\\p^\\lo\\Ti\\l < 

and 

\\kN,t\\l - ^ j dxdy\NujpN{x - y)\‘^ 

X [Pf (P + y)/2)Ppf (P + y)/2)P + P7 ((x + y)/2)7 

< c [|ptii^2|pf Hi + yf\\%i] < 
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The bounds (iii) follow from (see [6]) 

\\PN,th < \\kN,th cosh(||/cAr,t|| 2 ), and \\sN,t\\2 < l^N.th sinh(||/i:Ar,i||2) 


□ 


Lemma C.4. Let V be as in Lemma 2.1 and be defined as in (C.l), with ip = 
G (so that for all f G M, by Proposition B.l). Then 

(i) sup \\kN,t{-,x)\\2 

{a) sup \\pN,t{'^^)h, sup ||rAr,t(-,x)|| 2 , sup pAr,i(-,x )||2 < Ce^^^^ 

xGRS xSR® xGR^ 

(iii) sup \rN,t(.x,y)\ < Ce^^^^ 
x,y£K.^ 

Proof. The bounds can be proven as in Lemma C.l, taking into account the fact that 
the kernel of few,* is similar to the one of k^^t, just with p^{{x + y)f2y replaced by 
2 p^{ix + y)/2)ip^^{{x + y)/2). ’ □ 

Lemma C.5. Let V be as in Lemma 2.1 and k]\fp be defined as in (C.l), with p = 
p^Q G (so that for all t G M, by Proposition B.l). Then 

(i) ||Vpw,t|| 2 ,||Virw,P 2 <Ce^l'l, 

(m) ||AirAr,t||2 , ||AipAr,t||2 < 

Proof. Also in this case, the statements can be proven similarly as in Lemma C.2, keeping 
in mind that the kernels of pw,* aud r^v,* look exactly like the kernels of pw,* and r^v,* 
with the only difference that in the series expansion defining we have to replace one 
factor of A: 7 v,t or kj\f^t by its time-derivative (which simply means that in this one factor, 
we have to replace p^{{x + y)/ 2 )^ by 2 p^{{x -|- y)/2)p^{{x + y)/2)). □ 

Next, we establish the convergence of the kernel kM,t towards its limit, as A —)• oo. 
We denote 


h{x,y) = -a;f^”P(x - y) pi{{x + y)/2) 


with given by 


asymp 

UJ^ 


(x) 


bo 1 3 x^ 

^ R ~ ^ ^ 


for all X G with |x| < £ and by R^^™^(x) = 0 for |x| > £. Moreover, we dehne 
St = smh(kt), pt = coslfikt) - 1, n = sinh(A;t) - h. 

From (2.5), we find that 


^N,l{x) = 1 - 


sin(AR(PI - P) 


'N,e 


£cos(A]^_^(|x| - £)) 
\x\ 
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for all X € with RN ^ < |x| < i, where i? > 0 is the radius of the support of V. 
This gives the bound 

\NujnAx) - (C.8) 

for a constant C depending on i, and for all x S with |x| > RN~/^. 

We use (C.8) to bound the difference kN,t — h- 

Lemma C.6. Let V be as in Lemma 2.1 and kM,t be defined as in (C.l), with (p = 
e (so that for all t ^ R, by Proposition B.l). Let 

6 = min(/3/2,1 — /3). There exist constants C,ci,C 2 > 0 such that 

\\kN,t - hh, \\pN,t -Pth, \\rN,t - nh < CN~^eyLp{ci exp(c 2 |t|)) 
l|V {pN,t -Pt) lb, ||V(rAr,i - rbib < CN~^exjp{ci exp(c2|t|)) 

11^ iPN,t -Pt) lb, ll^(^Af,4 - n)lb < CA^“‘^exp(ci exp(c 2 |t|)) 


Proof. We estimate 

\\kN,t -kt\\l<2 j dxdy \Nujn,£A - v) - - y)|^ \Pt {A + y)/2)|" 


+ 2J dxd2/|a;f^™‘’(x-y)b {{x + y)/2)‘^ - pt{{x + y)/2)‘^ 

= I + II 


(C.9) 


Using (C. 8 ), we find 

|I| < I dxd2/^fc^^|(^f((x + y)/2)b 

+ C [ dxdy- — 2 |v?f((x + y)/2)b 

J\x—y\<RN ^ y\ 

< + CN-f^ 


On the other hand, using Proposition B.l, we find 


|II| ^ C* / dxdy 


/' 


l(|x-y| <i) 




<C dR 


ix - y |2 

pf^iR)^-ptiR? 


Tt{{x + y)/2f - Pt{{x + i/)/2)' 
2 


<C I dR\p^{R)-ptiR)\mp^{R)\^ + \pt{R)\^) 


< \\p^ - ptWl 

< exp(ci exp(c 2 |t|)) with 7 = min(/3,1 — (5) 


Hence, we conclude that 

II^AT,* - kt \\2 < exp(ci exp(c 2 b|)) 
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as claimed. Furthermore, we have 


PN,t -Pt = '^ '^{kthy {kN,tkN,t - ktkt){kN,tkN,tT~^~^ (C.IO) 

n>l ^ k j^Q 

This implies that 

\\PN,t -Pth < < CN~^exp{ci exp(c 2 |t|)) 

Similarly, we find 

\\rN,t -nh < C’iV“^exp(ci exp(c 2 |t|)) 

From (C.IO), we find 

V {pN,t - Pt) = ^{kN,tkN,t - kth) [ ^ ■7^-y{kN,tkN,tT~^ 

n>l^ 

- n—1 

+ v{hkt)[J2-7^^ J2^hkty-\kN,t~kN,t - kt~ktykN,tkN,tr~"~^ 

n>2 ^ ^k 

Therefore 

{pN,t - pt)h 

< \\VikN,t~kN,t-kth)h 

n>l ^ k 

^ 

+ \\V{ktkt)UkN,tkN,t - kkh E 7^ E 11 ^* 112 ' \\kN,t\\T-"-^^ 

n>l ^ ^k j=i 

< \\VikN,tkN,t - A:tfei)||2ell'=^.*ll2 + \\kN,tkN,t - hhh \\V{hkyh 

< C\\v{kN,tkN,t - ktki)\\2 + CN~^exp{ci exp(c2|t|)) 

We have 

{kN,tkN,t - ktkt)\\2 < \\'^kN,t{kN,t - h)\\2 + \\'^ikN,t - kt)kth (C-12) 

On the one hand, we find 

\\VkN,tCkN,t-h)\\l 

= j dxdydzidz 2 \S/ xkN,t{x, Zi)\ \{kN,t - kt){zi,y)\ \S/xkN,t{x, Z2)\ \{kN^t - kt){z2.,y)\ 

< j dxdydzidz2\{kN,t - kt)izi,y)y\VxkN,t{x,zi)\\VxkN,t{^,Z2)\ 

< C\\k]\f^t ~ kt\\2 

(C.13) 
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On the other hand, we compute 

- kt)h\\l= j dxdy J dzVx{kN^t{x,z) - kt{x,z))kt{z,y) 

< I dxdy / dz (Vx + yz){kN,t{x, z) - ktix, z))kt{z, y) 
+ I dxdy / dz {kN,tix, z) - kt{x, z))V^ktiz, y) 
integrating by parts. As in (C.13), we can estimate 


(C.14) 


dxdy 


dz {kN,t{x, z) - kt{x, z))V:,kt{z, y) 


< C\\kj^^t — /cflll 


As for the first term on the r.h.s. of (C.14), it can be estimated by 


/ 


dxdy 


dz 


Nu}n,i{x - z)ip^{{x + z)/2)Vip^{{x + z)/2) 

- - z)ipt{{x + z)/2)Vipt{{x + z)l2) kt{z,y) 

< sup [ dy\kt{z,y)\‘^ 

dxdz \Nu}]sf^e{x — z) — — z)f l^p^dx + z)/2)\^\S/ip^ {{x + z)/2)\^ 

+ J dxdz - z)f \<Pt{{x + z)/2) - (pt{{x + z)/2)^ ((x + z)/2)\^ 

+ J dxdz - z)f \(pt{{x + z)/ 2 )p|Vvjf ((x + z)/2) - Vipt{{x + z)l2)'\‘ 

Proceeding as in the analysis of (C.9) and using Proposition B.l, we conclude that 
/ 1/ dz{Vx + ^ z){kN,t{x,z)-kt{x,z))kt{z,y)^ <CN ^ exp(ci exp(c 2 |t|)) 

Inserting the last bound and (C.13) in (C.12), we find 

\\'^{kN,tkN,t - ktkt)h < C'iV“^exp(ciexp(c 2 |t|)) 

Hence, from (C.ll), we conclude that 

{PN,t - Pt )\\2 < C'iV“‘^exp(ci exp(c 2 |t|)) 

Analogously, we can show the bound for ||V(rAr^t — rt)||. 
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Also the bounds for ||A(p 7 v,t — Pt)\\ and ||A(rAr^j — rt)\\ can be proven in a similar 
way. In fact, observing that 

A^kNA^,y) = - naujnA^ - y){'Pt{{x + y)/‘^)f 

- 2NVu:n/{x - y)<Pt{{x + y)/2))Vipt {{x + y)/2) 

- Nujn,i{x - y){Vipt{{x + y)/2)f 

- NujnA^ - y)y^t{{x + y)/2)Ay5f ((x + y)/2) 

and that 

-nAun/A - y){y^t^{{x + y)/Af = ^n^^v{n^{x - y))fN,iA - y) 

- NXN,dN,£A - y)i(|a^ -y\ <i) 

we find 

sup f dy\AdN,td^y)\ 

X J 

and therefore, similarly to (C.13), 

\\AkN,t(kN,t -h )\\2 < Ce^''AkN,t - hh < C'A“^exp(ci exp(c 2 |t|)) (C.15) 


Moreover, integrating by parts twice, we obtain 

\\A{kN,t-h)h\\l< j dxdy j dz{V^+ Vzf{kN,tdiA - h{x,z))kt{z,y) 

+ I dxdy / dziVx + '^z)ikN,tix, z) - kt{x, z))Vzkt{z, y) 
+ I dxdy / dz{kN^t{x, z) - ktix, z))Azkt{z, y) 

This implies that 

||A(A:Ar,t - hyktWl < C'A"‘^exp(ciexp(c 2 |t|)) 
and thus, together with (C.15), that 

\\A{kN,tkN,t - ktk£)h < C'A“^exp(ci exp(c 2 |t|)) 


The last bound allows us to conclude that 

\\^{PN,t -Pt)\\2 < C'A“‘^exp(ci exp(c 2 |t|)) 

||A(rAry - rf )||2 < exp(ci exp(c 2 |t|)) 

□ 
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